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"EQUALIZATION METHOD AND EQUALIZER FOR CELLULAR MOBILE RADIO 
SYSTEMS " 



FIELD OF THE INVENTION 

The present invention relates to the field of the cellular telephone systems and in 
particular it relates to a space and time, quick convergence equalization method, for 
the cancellation of stationary and non stationary isofrequential interferes, in signals 
received from array antennas of base transceiver stations of a cellular mobile radio 
base system, and relevant systolic array equalizer. In particular, the invention refers to 
mobile systems employing a brief training sequence in transmission bursts. 

In telephone mobile systems, the main limitation to the traffic capacity offered to 
the public service is represented by the scarce band available. Therefore, all 
technological solutions enabling to increase the above mentioned capacity, or spectral 
efficiency, maintaining a given quality of the service in a wide area are particularly 
appreciated. To this purpose, appropriate cellular coverage strategies have been 
developed, together with multiplexing techniques in FDMA frequency (Frequency 
Division Multiple Access), and/or in time TDMA (Time Division Multiple Access), and 
more recently in space SDMA (Space Division Multiple Access) thanks to the use of 
the so-called intelligent array antenna combined to the receivers. 

The cellular coverage technique consists in dividing the service area in zones, 
called cells, each one served by its own base station BS (Base Station). Each cell 
avails of its own group of radio channels that can be reused in other duly spaced cells, 
in fact, the attenuation of radio signals is proportional to the distance with an exponent, 
or propagation coefficient, in general higher than 3. The whole of the cells employing 
the radio channels available to the system, without any reuse of the same, is called 
cluster of cells. Being the spectral efficiency inversely proportional to the cell area, 
microcellular systems are the most common approach to increase the spectral 
efficiency in high traffic density urban or suburban areas. As the cell area decreases, 
the problems associated to continuous handovers become more and more critical, 
consequently, mixed structures (microcells overlapped to umbrella macrocells) are 
generally preferred to simultaneously control quick mobile stations MS (Mobile Station) 
in scarcely populated areas, and slow mobile stations in high density areas. 
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In order to increase the spectral efficiency, a first approach as an alternative to 
the reduction of the cell dimension, consists in the reduction of the number of C cells-- 
per cluster, where C is also called cluster size, which is equivalent to increase the 
fraction of the total spectrum controlled by the single BS. This strategy is known as 
RCS (Reduced Cluster Size), whose theoretical limit is to render available the whole 
spectrum to all the cells, performing clusters of one cell only. 

A second approach consists in allocating a higher number of users in the 
spectrum already available to the cell, allocating more mobiles on the same channel, 
characterized by the same carrier and same time slot and therefore, without any 
intervention of the duster size. This second strategy is known as SCFR (Same Cell 
Frequency Reuse) and is based on the space division multiplexing SDMA within the 
cell. In practice, it is implemented in two different ways, the first one consisting in the 
division of the cell in smaller sectors through an arrangement of antennas with fixed 
radiation diagram (in general with 120° aperture), and a second and more advanced 
one, consisting in fitting the BS of relevant multibeam antennas, or synthetic aperture 
ones, that is capable to transmit and receive in the direction of mobile terminals with 
comparatively narrow directive beams; this implies that said beams must follow a 
generic mobile user moving in the azimuth plane. 

The formal difference between the use of an SCFR technique and the division of 
the cell itself in smaller sectors, resides in the capability to allocate the users in the 
relevant channels. Through sectorialization it is possible to manage more users on the 
same channel only if their position on the territory falls under different fixed sectors of 
radio coverage. In a scenario characterized by users in movement, this involves the 
need of continuous handovers. On the contrary, an SCFR structure is able to allocate 
more co-channel users in the same cell, without the need of partitioning the cell and 
have recourse to frequent handovers. Of course, even in this last case, the BS must be 
able to distinguish the co-channel users one from the other, in order to form 
appropriate radiation diagrams in receipt, each one capable to trace the relevant useful 
signal and reject the isofrequential interferes, corresponding to the remaining co- 
channel users simultaneously active within the same BS. and of course served by a 
same number of beamformers of. the same BTS (Base Transceiver Station). As we 
shall see below, the discrimination between useful and co-channel interferers in 
GSM/DCS systems takes place availing of a training sequence made of 26 known 
symbols forming the midamble of the transmission bursts, assuming to have assigned 
different training sequences scarcely correlated among them to different co-channel 
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users. The above mentioned sequence is purposely foreseen for the extraction of the 
symbol synchronism for the estimate of the channel impulse response. In the case of 
co-channel users of the same cell, they are also synchronized among them, and 
therefore it is possible to locate the same through analysis of the training sequence 
only. 

When adopting the RCS and SCFR strategies, it is absolutely necessary to avoid 
that co-channel interferers can worsen the quality of the demodulated signal, thus 
making vain the advantage of the higher capacity of enabled traffic. We want also 
underline that, in the systems employing an approach of the RCS type, the interference 
» is due to the co-channel users of adjacent cells, while in systems employing the SCRF 
approach the interference is generated both by the co-channel users of the same cell, 
and by the co-channel users of adjacent cells. 

The mobile radio channel has been studied for a long time, and different models 
have been proposed, consistently with the more frequent transmission environments, 
5 rural, urban, etc. From these studies it resulted that, due to the space and time 
dispersive nature of the radio channel, in addition to the interferers generated by co- 
channel users on the direct path, there are also those corresponding to the multiple 
paths. 

Particular attention shall be paid to possible delayed echoes of the useful signal. 
0 since in some examples of the known art they are confined to the role of co-channel 
interferers, while in the embodiment according to the present invention the same are 
actually recognized as integral part of the useful signal. 

Our attention shall be concentrated below, on base stations employing antenna 
arrays to implement the RCS or SCFR techniques, the latter also according to the more 
25 advanced SDMA technique, with dynamic assignment of the channels, or concurrently 
with the use of traditional techniques of fixed sharing of frequencies. 

We shall also start from the assumption that the synchronization problem has 
already been solved, that is the acquisition method of the carrier synchronism is known 
as well as the position of the training sequence inside the burst. 
30 BACKGROUND ART 

The nature of the mobile channel is such to cause a spatial dispersion (more 
arrival directions of echoes) and time dispersion (delay of echoes) of the power of the 
signal reaching each antenna. In a context characterized by propagation along multiple 
paths, it can be reasonably expected that the different paths of the signal reach an 
35 antenna array with different DOA. This characteristic of the mobile channel can be 
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used to perform an equalization in the spatial dimension, it is possible for instance to^ 
filter only the maximum power path, or appropriately combine the signals impinging 
from different paths among them to cancel co-channel interferers. The equalization of 
time dispersion shall be considered hereafter. 

5 In order to be able to effectively discriminate the useful from the different kinds of 

interfering signals, the array must be able to dynamically modify its own directivity 
function and insert some zeroes in its own radiation diagram in the directions of the 
interfering signals. We shall give below the description of some interfering implemented 
according to the known art, highlighting time by time the relevant drawbacks. The 

10 straightforward description, which however nothing detracts to the general character of 
the same, is limited to a linear and uniform array of antennas, called also ULA (Uniform 
Linear Array). 

Making reference to fig. 1, the array consists of M equal aligned and equispaced 
elements. The separation between adjacent elements is d. For narrow band signals 
15 with frequency of the carrier fc, the corresponding wavelength can be defined as: 

A = 7- (D 

J c 

where c is the propagation speed (c = 3.108 m/s in vacuum). The spacing d between 
the elements can be conveniently expressed according to the wavelength A as: 

rf = AxA (2) 
20 A is thus the separation between the elements in terms of wavelength. 

Assume later on, that a sole punctiform source is present, sufficiently distant to 
render the wave impinging on the array similar to a plane one. The direction indicated 
in fig. 1 by angle 0 is the direction of arrival of the plane wave, hereinafter called DOA 
(Direction Of Arrival). Therefore we obtain that the relative delay r between the signals 
25 received by two adjacent elements of the array is: 

doos{0) 



x = • 



(3) 



c 

the real signal received on the k-th element of the array is then: 

r k (t) = s(t~T k ) = s(t-kr) (4) 
It is assumed that the delays relevant to the different receivers along the array 
30 are small compared to the reciprocal of the signal band, so that the same can be 
approximated with opportune phase variations of the complex envelope r k (t) . 
Therefore it results: 
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r k {t) = s{t)e-W (5 ) 
the phase factor2;r/;r can be written in different ways: 

l7tf t = 2rf c = 2tt^co^0) = 2;rA cos(0) (6) 

Finally, if ?(/) is the equivalent in base band of the signal received on the first 
element of the array, and defined a vector a(0) as: 



■(*)« 



i 

- jlitf*. cos(0) 



e -A*{ M-l)Acos(0) 



(7) 



we can group the signals received on the different elements of the array in a vector 
r(r): 

?(/) * *(e)7{t) 

Vector a(0) is the so-called steering vector, or array response vector, at an 

impinging plane wave with DOA equal to 9 . 

In the background art, one of the first approaches to the problem of spatially 
selective reception and transmission is that of the so-called Switched Beams systems, 
that is those employing antenna arrays to form a fixed number of beams oriented in 
fixed directions. Beams so obtained, are selected, both in transmission and in 
reception, in order to employ the better oriented beam in the direction of the mobile 
terminal. Though easy to implement, these systems have the great disadvantage not to 
be able to effectively discriminate the interfering signals. 

Making reference to fig.2, we notice a Switched Beams receiver including a four- 
element array, connected to the inputs of a Beamformer block, whose four outputs are 
connected, to the inputs of a beam selection block and to those of a selector, 
respectively. The output of the beam selection block is connected to a control input of 
the selector. 

During the operation, the beamformer of fig. 2 is used to form a radiation diagram 
consisting of a fixed number B of beams in pre-set directions & on the azimuth plane. 
The inputs of the beamformer are represented by the vector r(f) (hereinafter denoted 
with y(r) ), neglecting for convenience, the tilde sign of the complex envelope. 

The beamformer consists of a matrix of coefficients [w mfA ] defined as follows: 




with m = 1 , 2 M and b = 1 , 2 B. The signal on the b-th output of the beamformer 

shall therefore be equal to: 

y ^ (,) = ^ 5 ( r ) e -y2»(«-0A(cni(^)-cos(fl)) ^qj 

If 0 = ^ the output results the maximum one. In case more signals are present, 
impinging from several directions, considering the arrival direction # = the 
contributions of the different elements of the array are summed up in phase. The b-th 
output of the beamformer shall therefore be characterized by a radiation diagram with 
maximum in the direction ^ . 

The selection of the best beam takes place according to a criterion tending to 
minimize the error probability, that is applying a MLSE estimate method (Maximum 
Likelihood Sequence Estimation) of the sequence transmitted at maximum likelihood . 
The above mentioned estimate is profitably made applying the well known Viterbi 
algorithm to a trellis sequential diagram, or trellis, whose branches represent all 
possible status transitions of the receiver. The algorithm calculates, for all the possible 

sequences {x„} , the distance between the corresponding signal transmitted: 

^lk}) = Z„^-«r) (H) 

and the received signal y{t) , and then selects the sequence {x n } corresponding to 
the minimum one among said distances. From the above we can understand that the 
higher are the distances among the possible signals >{f| {*„})■ the lesser is the 
probability to select the wrong sequence during the demodulation. Therefore, defining 
^min the minimum value of the distance among all possible y{t\ {*„}), this parameter 

results being a reliable index of the error probability. 

Due to the calculation complexity of this approach, in some known examples 
another method is also proposed for the beam selection, simply based on the quality of 
the demodulation of the training sequence. Among all the possible beams, it is selected 
in fact the one determining the smallest number of errors in the midamble 
demodulation. At equal number of errors made, it is used the beam that, at the end of 
the Viterbi algorithm applied to the midamble, originates the minimum distance. 

An approach giving the antenna a higher intelligence level versus Switched 
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Beams methods, essentially consists in the execution of a spatial filtering on the 
signals received by the antennas, by a beamformer resuming in the output signal the 
directivity characteristics of an opportune radiation diagram capable to reduce the co- 
channel interference. This approach, in the major part of the known application 

5 examples, required the prior knowledge of the DOAs of the useful and co-channel 
interferers. Different types of algorithms studied to solve the problem of the DOA 
estimate are found in literature, availing of an array of receiving antennas. The most 
representative algorithm is that known as MUSIC (Multiple Signal Classification), 
proposed by R. O. Schmidt in the paper under the title "Multiple Emitter Location and 

10 Signal Parameter Estimation", published in IEEE Trans. Antennas and Propagation, 
vol. AP-34, No. 3, March 1986, reference to which shall be made for a more complete 
description. The algorithm requires the restrictive assumption that the covariance 
matrix of the signal is not singular; in this case the signals on antennas are not 
correlated among them and the steering vectors result linearly independent. 

15 When arrival directions are known {§ n / = 1, ,p) an appropriate radiation 

diagram shall be constructed. A possibility is to define as many diagrams as are the 
estimated DOAs (p) f on the basis of the MVDR technique (Minimum Variance 
Distortionless Response), described for instance by S. Haykin in the volume "Adaptive 
Filter Theory", published by Prentice-Hall, Inc., Englewood Cliffs, N.J., 1991. 
20 Said approach enables to estimate the beamformer coefficients in order that the 

radiation diagram has unitary response coinciding with the selected direction, 
minimizing at the same time the interference effect and of the noise coming from all the 
other directions. 

Therefore, up to this stage the arrival directions of the different signals have 
25 been estimated, making for each one of them an appropriate radiation diagram. 
However, according to the described technique it is not possible to know which is the 
arrival direction associated to the signal considered as useful and consequently, which 
is the radiation diagram more convenient to be used to perform the demodulation. One 
can think to select the beam according to the methods previously described, in order to 
30 select the one that enables a better demodulation of the training sequence. 

In the case of DOA estimate, the above arguments are strictly valid for the case 
of p signals uncorrected with p<M, where M is the number of elements of the array, or 
for scenarios characterized by a low number of multiple paths, as the mobile 
communication systems in rural area. 
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In the case of propagation in environments characterized by number of multiple^ 
paths, such as for instance in metropolitan areas, the extraction of the DOAs can result 
ineffective, since the number of arrival directions of the different signals is often higher 
than the array size, and also the hypothesis of uncorrelation among the different 
5 signals on the antennas is scarcely verified. However, at least in the case where 
multipaths relevant to a same signal belong all to an angular sector, not too wide, the 
performances of the algorithm often remain satisfactory, and in any case higher than 
those guaranteed by the Switched Beams approach. 

To render the direction estimate more efficient in case of correlated signals, one 

10 can think to use other DOA research algorithms, like the WSF (Weighted Subspace 
Fitting) proposed by B. Ottersen and M. Viberg in the article "Analysis of subspace 
fitting based methods for sensor array processing", following the Proc. ICASSP 89, pp. 
2807-2810, Glasgow, Scotland, May 1989. This algorithm, at the cost of a higher 
computation complexity, is able to process also coherent sources. The advantages are 

15 however highly limited as already said both by the large number of arrival directions 
and by the fact that, even succeeding in estimating more than one direction for each 
signal, it should be necessary to use an algorithm for the beam composition capable to 
completely exploit the knowledge of said directions. Availing of an effective 
demodulation diagram, capable to control the multipaths, as the one based on an 

20 MLSE algorithm, it can result useful to point the beam toward more replicas of the 
same signal, in order to maximize the signal-noise ratio, thing which is not possible with 
the already mentioned MVDR technique which, as already said, considers one sole 
direction as useful. 

The algorithms described above, foreseeing the DOA estimate prove also to be 
25 ineffective in case of non-stationarity during transmission. Let's consider in fact the 
case of an interferer disturbing the transmission only in a semiburst (and not during the 
midamble). The arrival direction of an interferer, estimated using a beam selection 
algorithm based on the quality of the midamble demodulation, could not consider the 
effect of said interferer. In this case, it would be necessary to have recourse to more 
30 complex algorithms for the selection, or composition, of the beam based on the 
demodulation not only of the midamble but of the whole burst. 

In the cellular mobile system, the number of interferer signals may sometimes 
result higher than the number of the antennas forming the array, in this case the 
hypothesis p < M is no more verified and the array can also not be able to insert zeroes 
35 in the radiation diagram in the direction of each interferer. This means that the methods 
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based only on the DOA estimate additionally loose their effectiveness. However, ^ 
also in this context, the adaptive array could be able to combine the signals presents at 
the antennas to improve the ratio between the signal and the sum of interferers with the 
noise, ratio known in literature under the acronym SINR, (Signal to Noise and 
5 Interferers Ratio). In the following articles by J. H. Winters: 

♦ "Optimum Combining in Digital Mobile Radio with Cochannel Interference", 
published in IEEE J. Select. Areas Commun., vol. SAC-2, No. 4, July 1984; 

• "Optimum Combining for Indoor Radio Systems with Multiple Users", published in 
IEEE Trans. Commun., vol. COM-35, No. 11, November 1987, is developed the 

10 concept of optimum combiner. 

This is an alternative technique to the DOA estimate, which bases on the "prior" 

knowledge by the receiver of a sequence of known symbols {*„} (training sequence), 

present in the single transmission bursts for synchronization purposes. The sequence 

[x n } is additionally exploited in the spatial filtering process, through which the signals 

15 received from the antenna array are appropriately combined to maximize the SINR. 

This approach enables to avail of the spatial difference of the signal, besides 
reducing the relative power of the co-channel interferers, also to face the Rayleigh 
fading of the useful signal. 

Fig. 3 is the block diagram of an M elements optimum combiner. As it can be 

20 noticed, the combiner has M inputs for a same number of base band signals 
^W.-JiWr-JAfW' received by corresponding M elements of an array (not 
shown in the figure). The above mentioned signals are multiplied by the relevant 
weights w t ,... 9 w i9 „. 9 w M , in general complex. Weighted signals reach the inputs of a 
first adder, at the output of which a y out i n ) signal is present, coming out from the 

25 combiner and reaches also, with negated sign, an input of a second adder, at the other 
input of which a x n signal is present. At the output of the second one, an error signal 
e n is present, reaching a WEIGHT ADJUSTMENT block having other inputs for signals 
y y (n) 9 .„ 9 y { (n) 9 ... 9 y M {n). At the outputs of said block, weights w l9 ... f w n ... 9 w M are 
present. 

30 In the operation, the signal y^M can be expressed as: 

y 0 M = ty t (n)w t (12) 

jo I 
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Denoting with w = [m> m w 2 ,...., w m ] T the vector of the weights and with 

y(n) = [y l (n), y 2 {n),....,y M {n)] T the vector of the signals received and adopting a 
matrix notation, we can define the weighted sum of the signals coming from the 
antennas as: 

^W = w T y(«) (13) 
The difference between the desired signal x n , corresponding to the samples of 
the training sequence, and the signal y out (n) at the output of the combiner defines the 
error signal: 

This signal is used as control signal in the adaptation process of the combiner 
weights. 

The scope of the adaptation process is to find a vector of weights w, such to give 
at output, at any instant, a signal equal to the desired one x n . When the input signals 
can be seen as stationary stochastic variables, the WEIGHT ADJUSTMENT block 
searches the vector of the weights minimizing a cost function corresponding to the 
mean quadratic error MSE (Mean Square Error), defined as follows: 

«'(w) = E[\e n \ 2 ] 05) 

where £[•••] indicates the average value of the subject in brackets. 

By replacing the (14) in the (15) we can obtain the whole of the weights 
minimizing the MSE: 

e 2 (w) = e\\x J' +w*y(«)y"(n)w-2w'x n y(«)] 

I- (16) 

= £[I*J 2 ] + W ' S » w - 2w ' 5 *c 
where the spatial covariance matrix R w of signals at the antenna input can be 
expressed as follows: 

R„ = E[y(n)y(n)} (17) 
and the cross-covariance matrix among the M signals at input of the antennas 
and the desired signal, can be expressed as: 

R„ = e[y(«K] < 18 > 
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The mean square error defined in the (16) is a quadratic function of weights w. 
The components of the gradient of function MSE, are the partial derivatives of the MSE 
versus the value of weights. Deriving the (16) versus w ( we obtain the gradient as 
linear function of weights: 

V£[|eJ 2 ] = 2R„w-2R^ (19) 
When the selection of weights has been optimized, the gradient results null, then: 

R " W °" = (20) 

The weight vector w opt is therefore the one giving the minimum mean quadratic 

error. The equation (20) is called Wiener-Hopf equation. For a more detailed analysis 
of the research methods of the optimum weight vector, reference shall be made to the 
article "Adaptive Antenna Systems B , by: B. Widrow, P. E. Mantey, L. J. Griffiths and 
BB. B. Goode, published in IEEE Proceedings, vol. 55, No. 12, December 1967. 

Considering now the time equalization, it is opportune to introduce the following 

considerations: let {x n } a sequence of N symbols transmitted on a channel with 
dispersion length L, that is whose impulse response is significant within a time LT, 
being T the symbol time. Let's consider an array of antennas consisting of M elements. 

The vector of signals received on the array of antennas y{t) =[y x (t), y 2 (t),:..-,y M {t)] 
is represented by the convolution of a vector h(r) with the transmitted sequence {x n } : 

yW=£* n h(r-*r) + n(r) (21) 

where the vector elements h{t - nT) are the samples of the impulse response of the 
channel on the different antennas. The vector h(r) can be written as follows: 

hto = X4e- /2 * r A('-r,MS) ( 22 ) 

being: 

• N the number of multiple paths, 

• A t the path loss of the i-th path, 

• 0. the arrival direction, 

• h Q the convolution of the modulation impulse and impulse response of the reception 
filter, 
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• a(6J.) is the steering vector, or vector of the array response in direction 0 {% defined 
in the (7). 

The vector n(r) represents the additional noise, we shall time and spatially 
consider white and that satisfies the following conditions: 

5 r J (23) 

E[n k {t)n k {t + r)] - a 2 S(r) 

In the major part of the known realization examples of an optimum equalizer for 
the compensation of distortions introduced by the channel on time characteristics of the 
sequence transmitted, the equalizer is made up of a filter adapted to the channel 
impulse response, followed by a maximum likelihood detector (MLSE) based on the 

10 Viterbi algorithm. Both for the realization of the adapted filter, and for the MLSE 
detector one, it is absolutely necessary to known the channel impulse response, 
generally unknown. Availing of the knowledge of the midamble, the equalizer estimates 
the channel impulse response, to calculate then the coefficients of the filter matched to 
the channel impulse response and be able to estimate the sequence transmitted, 

15 profitably using the Viterbi algorithm. The possibility to perform an estimate of the 
vector h(/) makes therefore useless the explicit calculation of the same according to 
the (22), which presupposes the estimate of the arrival directions. 

Algorithms mentioned up to now, face the problem of the estimate of the 
sequence transmitted in two different equalization phases, namely: 

20 • spatial equalization for the cancellation of interferers; 

• time equalization of the spatially filtered signal. 

A separate equalization enables to spatially combine only the echoes of the 
signal resulting not significantly delayed. The delayed echoes of the signal are in fact 
considered as co-channel interferers. However, this approach can result in a reduction 
25 of the ratio between signal and noise at the time equalizer input, as it shall be noticed 
through comparison of figures 4 and 7, and in a lower robustness of the beamforming 
algorithm versus particular azimuth distributions of echoes, as it could be noticed 
comparing figures 5 and 8. 

In the light of what just said, let's consider the cost function (15) of the 
30 Wiener-Hopf combiner, previously dealt with: 

e 2 (w l ,w 2 ,....,w M ) = E^e„\ 2 ] = E 



(15) 
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where w, is the coefficient according to which the signal received at the i-th antenna is 
weighed. This algorithm maximizes the ratio between the signal and the noise, where 
signal means all those echoes resulting correlated with the training sequence {x„ } . If 
the channel introduces delays higher than the symbol time, T = 3J/js for GSM, 
5 echoes so delayed result decorrelated with the training sequence and are treated by 
the combiner as co-channel interferers. 

Fig. 4 is the radiation diagram on the azimuth plane of an antenna array, where 
the beamforming coefficients are calculated with the Wiener-Hopf algorithm applied to 
the functional (15), in presence of an useful signal generated by a mobile MSd, with 

10 relevant echo delayed by 5 jxs, and of a co-channel interferer generated by a mobile 
MSi. As it can be noticed, the radiation diagram shows a maximum in the useful 
direction (90°), a first zero in the co-channel interferer direction (120°), and a second 
zero in the useful signal echo direction (60°). Due to the second zero, the possibility to 
avail also of the power of the useful signal in the estimate process of the transmitted 

15 sequence is lost. 

Fig. 5 is the radiation diagram on the azimuth plane of the antenna array used in 
a scenario differing from the previous one only for the presence of a replica of the co- 
channel interferer placed in the same arrival direction (90°) of the main signal. As it 
can be noticed, the radiation diagram still shows the two zeros in the previous 

20 directions, of the interferer (120°) and of the echo of the useful signal (60°), 
respectively, but in this new scenario, the maximum in the useful direction (90°) is no 
more present. This behaviour can be explained, considering a traditional 
synchronization system, by the fact that the combiner sees a unique useful signal, 
different from the interferer, in the common direction (90°), however attempting to 

25 reduce the effect of the interferer in that direction reducing the gain. This involves an 
additional reduction of the power of the useful signal present in the estimate process of 
the transmitted sequence. Therefore the combiner looses the opportunity to transfer 
the maximum gain lobe in the sole echo direction of the useful signal (60°), thus 
succeeding to nullify the echo effect of the interferer having the same direction of the 

30 useful signal. 

In mobile radio environment more recent algorithms than those described up to 
now are known, that try to overcome the drawbacks highlighted in the description of 
figures 4 and 5. The above mentioned algorithms still avail of the concept of 
minimization of the mean square error expressed by an appropriate functional, like it 
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occurred for the optimum combiner, but, contrarily to this one, the information on the 
channel included in the expression of the above mentioned functional is now given by a 
so-called spatial-time filter due to the fact to adapt to the total response of the channel. 
A first example of the above is represented by the international patent application 
5 WO 97 08849 A, published on 6/3/97 under the name of THOMSON CSF, under the 
title: "METHOD AND DEVICE FOR THE SPATIAL MULTIPLEXING- 
DEMULTIPLEXING OF RADIOELECTRICAL SIGNALS FOR SDMA MOBILE 
SYSTEMS", with French priority of 22-8-95. 

In the above mentioned application it is described a base station employing in an 

10 "intelligent" manner an array of antennas to which a multiway transceiver is coupled, 
equipped with beamforming means particularly suitable to enable the assignment of a 
same frequency-time channel to more mobiles in the same cell, provided that the 
relevant signals have sufficiently decorrelated arrival directions. The invention, which 
essentially claims an application of the SDMA technique in mobile radio environment, 

if shows two embodiments distinguished by different typology of the beamforming means 
used in reception. More in particular, the first embodiment is referred to means 
requiring a functional block for the estimate of the steering vectors, which are then 
passed to a true beamforming spatial filter. The beamforming means of the first 
embodiment can be compared to those already mentioned for the DOA estimate and 

20 result far less adequate to the present invention compared to the beamforming means 
foreseen in the second embodiment, therefore the first embodiment shall not be 
examined. However, it is useful to remember that the techniques used in the first 
embodiment are called "transparent" because they do not require any prior knowledge 
of the sequence of symbols received. 

25 Beamforming means foreseen in the second embodiment of THOMSON paper 

use, on the contrary, the prior knowledge of the training sequence, as it occurs in the 
invention described below. These second means include a space-time filter. The 
optimization of the total response of the channel is made through two algorithms 
alternative between them, called FAST-BBT-DFE and FAST-DFE, respectively, where: 

30 • FAST is an acronym of Space-Time Matched Filter, 

• BBT is an acronym of Bruit Blanc Temporellement, with French terms, and 

• DFE of Decision Feedback Equalizer, with English terms. 

We quote below what is said on the advantages of these algorithms: 

• "The FAST-BBT-DFE enables to rephase the power of all the paths of the useful 
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signal and to reject in sub-optimal way the signals interfering from multiple paths. 

This procedure involves a gain in terms of signal to noise ratio on the current symbol 

compared to an only-spatial filtering algorithm (FAST-DFE), in fact the algorithm 

shows the same rejection capabilities of interferer signals with the difference that the 
5 echoes of the useful different from the main path are no more assimilated to 

interferer signals, so that a same number of freedom degrees are cleared to 

consider other interferers. 
# The FAST-DFE enables to process a signal with interferers coming from an array in 

an almost optimal way, since it performs not only the phasing of the echoes of the 
10 useful signal, but also that of interferers before their rejection. This procedure 

enables to eliminate an interferer source using one sole spatial freedom degree 

whichever is the number of echoes according to which this source is received by the 

array". 

The main restrictive hypothesis of FAST-BBT-DFE algorithm consists in 
15 considering the matrix K b as having impulse nature, which involves interferers 
completely uncorrected between them, and therefore comparable to a white noise in 
the band of the useful signal. In other words, the interferers are considered free of 
echoes, or better, each echo is considered itself an autonomous interferer. Under 
these assumptions, the convergence of the algorithm appears rather quick. 
20 It is well clear that the application of such an algorithm is not recommended in 

mobile radio environment where, particularly in urban small size cells, the typical 
scenario is on the contrary characterized by co-channel interferers fit with a range of 
relevant echoes from multiple paths. The fact not to consider this reality involves the 
disadvantage not to be able to optimize in an actually credible way the total 
25 signal-noise ratio, or, whenever the SCFR technique (in SOMA mode) is used in order 
not to be able to separate the different co-channel mobile users in an optimum way on 
the azimuth plane. 

The FAST-DFE algorithm, shown making reference to fig. 12 of the quotation, 
does not place any limit on the noise nature, including the co-channel interferers, thus 
30 treating the channel in an almost optimal manner. Also in this case the structure 
consists of a space-time filter associated to a decision organ of the DFE type, even if 
the use of the Viterbi is considered possible, without describing methods for said use. 
The space-time filter consists of N time filters W*, N HR long, one for each w-th way 
branching from the N antennas. An additional time filter HR is also present, it too N ^ 
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long, placed on the feedback path of the DFE equalizer. As it is known from the theory 
of equalizers, a DFE, before submitting the symbols received to a threshold decider, 
tries to remove the ISI from the sequence of the relevant samples, providing to this 
purpose a deleter filter, acting in feedback on postcursor samples of the impulse 
5 response, and of a preliminary equalizer acting forward on precursor samples. In the 
diagram of fig. 12 the HR filter represents the deleter, while the whole of the filters 
whose filtered signals are summed up to form a unique sum signal, can be compared 
to the preliminary equalizer. The signal cleaned by ISI, denoted in the text with z{n) f is 
obtained from the difference between said sum signal and the signal filtered by HR. 

10 The z{n) signal crosses the threshold decider that gives an estimate z{n) of the data 
sequence transmitted. The high number of coefficients to update involves a given 
slowness in the convergence of the estimate and therefore the need of long training 
sequences, this fact being autonomously highlighted in a table of the advantages and 
drawbacks of the FAST-DFE algorithm. 

15 If one wants to employ the above mentioned space-time equalizer in GSM 

environment, where the training sequences foreseen by the standard strictly have a 
lengths of 26 symbols, and can therefore be considered short, the problem arises 
about the convergence of the estimate within the time slot of such a sequence. In this 

case, in fact, the estimate of error £ would undergo an additional slackening due to 
20 lack of accuracy in the evaluation of the error itself, the DFE having to this purpose use 
the £ (n) estimates in place of the no more available samples of the training sequence. 
Therefore, there is the actual possibility that a great part of symbols of the whole 
transmission burst is wrongly estimated, causing by this a degradation of the quality of 
the demodulated signal. Also, in case of transmission systems adopting short bursts, 
25 which could also include the GSM itself, with a burst of 148 symbols only, there is also 
the risk of a missed convergence within the whole duration of the same burst. 

The inconvenience highlighted is in itself sufficient to discourage the use of such 
an algorithm in GSM environment, but additional considerations shall be made to point 
out the consequent technical problems of non negligible importance. 
30 A second problem that can originate from the slowness in convergence of the 

estimate is that of the poor capability of the equalizer to face the channel non 
stationariness, mainly due to interferers that enable outside the time slot of the training 
sequence, and/or to the case in Which the channel can no more be considered 
time-invariant, as it occurs for instance when the Doppler effect is significant, due to 
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the high speed of the mobiles. The above mentioned inconveniences become even 
more important in presence of sudden and high intensity interferers (c / / <-10 dB). 

This is a condition that can frequently occur in urban mobile radio environment, for 

instance, when due to a so-called "comer attenuation" the power of the interferers 
5 suddenly surpasses that of the useful signal. Particularly in this last context, more 

powerful solutions are required than that indicated in the (26), consisting in weighting 

each iteration of the estimate with the forgetting factor X < 1 . 

A last drawback, originated by the same causes of the previous ones, consists in 

the difficulty to physically realize a receiver able to perform in one sole symbol time all 
10 the arithmetic operations necessary to the adaptation of a so great number of 

coefficients of the space-time filtering structure. 

A second example in which the space and time term "jointly" enter the LMS 

algorithm is disclosed in the German patent application DE-A-1 963941 4, inventor: 

RADEMACHER L.. The claim 1 of the citation recites textually: 
15 Method for the parametrization of a receiving station (BS) comprising adaptive 

antenna devices (AE) for receiving received signals (es) via time-variant channels, in 

which method 

- the received signals (es) from at least two adaptive antenna devices (AE) are 
weighted with antenna coefficients (a) and are combined to form received data (z), 

20 - error magnitudes (e) are determined in a device for determining error magnitudes 
(e) from the comparison of received data (z) and reference data (y) modelling these 
received data (z), 

- channel coefficients (RW) being provided for modelling the reference data (y), 

- these error magnitudes (e) are supplied to a channel estimator (ACE), 

25 - the channel estimator (ACE) jointly determines channel coefficients (RW) and 
antenna coefficients (a) minimizing the error magnitudes (e), and 

- the channel coefficients (RW) are provided for a detection of the received data (z) 
to be evaluated. 

The detector is a Viterbi estimator whose short-time decisions are used outside 
30 the training sequence of a GSM cellular system data block. Each channel coefficient 
RW is the product of the effective channel coefficient W by a phase rotator term R, 
common to a number of channel coefficients W, foreseen to equalise a time-variant 
channel by compensating the Doppler shift due to the velocity of the mobile. The 
optimization LMS algorithm has recourse to three separate gradients in vectorial form: 
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a first gradient is for updating the weights a t , a second gradient is for updating the 
channel coefficient W t? and a third gradient is for updating the phase rotator term R t . In 
each gradient the updating term at the instant t-1 is multiplied by an its own scalar 
leakage factor L. The three leakage factors L a , L w , L R include a respective product <xn a , 
5 pnw. p|iR f where a, (5, p are constants established according to the principle of the 
steepest rise of the gradient, and nw, hr are adjustable adaptation step lengths. The 
effect of various combinations of \x at ^iw. hr is described in the text. It can be appreciate 
that all the entities, except error e, are separate groups of three. 

In the above method for parametrization of a receiving station, no solution in 

10 closed form is provided, but only the recursive strategy to move toward the optimum 
value of the parameters. The optimization is carried out with respect to each gradient of 
the objective function, one independently from the others, by using the same error 
(probably this originates the term "joint" ); consequently inter-relationships between 
spatial and temporal terms are neglected and the opportunity to exploit these inter- 

15 relationships to improve the equalization definitively lost. No proof of convergence is 
given, in any case the gradient method for solving LMS problems is universally known 
to have a slow convergence. In fact to track rapid variation the phase rotator R is 
provided. Besides a right calibration of the adjustable adaptation step lengths n, hr 
is necessary. A reasonably hypothesis is that the phase rotator R be not in grade to 

20 fight a strong non-stationary interferents, that because the effect of such interferents is 
not comparable to the continuous phase rotation of the Doppler shift, in case of the 
strong interferent the convergence wouldn't be reached in a short/medium time. 

PURPOSES OF THE INVENTION 
Main scope of the present invention is to overcome the above mentioned 

25 drawbacks having as main scope that to indicate a spatial-time equalization procedure 
particularly quick in converging towards optimum performances, even outside the 
midamble, and that can be profitably employed in base stations of a cellular mobile 
system, for instance GSM, fit with so-called smart antennas. The main scope of the 
invention cover a non-stationary channel, particularly if it is affected by sudden and 

30 high intensity interferers. assuring the shortest possible duration of the transient within 
which a loss in quality of the demodulated signal occurs and the limitation of such a 
loss to limited values. 

A second scope is to indicate an optimization algorithm whose solutions are 
obtained in closed form by Direct Matrix Inversion. 
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A third scope is to indicate an iterative digital method for the updating of an 
estimate of the optimum equalization parameters, accurate and of quick execution at 
the same time. 

Fourth scope of the invention is to indicate the realization methods through 
5 systolic arrays of an equalizer based on the just mentioned iterative digital method. 

SUMMARY OF THE INVENTION 
a) To attain these aims, subject of the present invention is an equalization method 
to cancel the isofrequential interferers in signals received by array antennas of base 
stations of a mobile radio system, as described in the first claim. 
10 Profitably the decider component is a Viterbi processor exploiting said optimized 

samples of the impulse response of the channel to calculate, in the known way, branch 
metrics of a trellis, used to estimate at maximum likelihood the sequence of symbols 
transmitted. 

In the method of the invention the combined estimate of the best equalization 

15 parameters exhausts at the last sample of the training sequence, generally included in 
a midamble of the transmission burst The 26 samples of the GSM midamble (out of 
the 148 available) are sufficient to assure the convergence of the estimate even in 
presence of medium quality signals received. The best equalization parameters should 
remain then available for the demodulation of the remaining part of the stored burst, 

20 and the demodulated signal does not loose quality, provided that channels are 
time-invariant and the interferers have stationary characteristics, but these ideal 
condition are not met in practical cases. 

The method according to the invention solves in the adaptive way the problem of 
a correct equalization in presence of channels with time-variant characteristics and/or 

25 non-stationary interferers, updating the above mentioned estimate of the best 
equalization parameters at each new symbol received. From the operational point of 
view, it is necessary to this purpose to avail both of a vector of antenna input signals, 
and of a vector of data transmitted. While the first vector is always available, the 
second vector is completely known only during the midamble, where it corresponds to 

30 the training sequence stored in the receiver. We see therefore the need, once the 
midamble is exhausted, to replace the training sequence by the symbols of the 
estimated sequence, considering that the symbols already decided at an instant n t are 
but those associated to a path selected in the trellis up to instant n-th. Simultaneously, 
an updated sequence of optimum samples of the impulse response of the channel is 

35 sent to the Viterbi processor, for the updating of the values of the branch metrics. This 
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way of proceeding inevitably introduces an excessive delay in the estimate of the 
optimum equalization parameters, due to the reaction time of the Viterbi algorithm 
which adds to the delay due to the time matched filter, thus preventing to follow quick 
variations of the channel characteristics. 

5 The inconvenience is avoided forcing the Viterbi algorithm to decide just after few 

steps (typically 4, or 5) which is the most possible sequence of data, that is taken at the 
time of the present symbol climbing in the trellis the maximum likelihood path, while 
final decisions for the MLSE detector are taken at the end of the burst, in order that the 
same are as much reliable as possible. 

10 Since short term decisions can result wrong, an error propagation phenomenon 

can be generated, which can last until wrong data are reused for the estimate of the 
best equalization parameters, that is for a number of symbol times equal to the intrinsic 
delay the short term decisions are made available. 

The equalization method of the invention, eliminates also this second drawback 

15 copying as a whole, at any symbol time, the L symbols necessary to the updating of 
the equalization parameters. It is well evident that the updating of the best equalization 
parameters is a new combined estimate of the same. 

Profitably, it is possible to proceed to the above mentioned updating, and to the 
consequent adaptive equalization, starting from the ends of the midamble and 

20 proceeding in parallel in the two directions towards the ends of the burst previously 
stored and synchronized. This involves the advantage to avail of the double of the time 
for calculations. It is also possible to proceed with a first pass on the right semiburst 
and on the left one, only afterwards, using the decided data of the right semiburst to 
have a longer group of data on the left one. This procedure improves the accuracy but 

25 not the speed of the estimate. 

With the actions made up to now, and considering that the combined estimate of 
the optimum equalization parameters is extended to the whole burst, it is useful to 
increase the effectiveness of the more recent data versus older ones, in order to render 
the method more reactive in respect with the variations of the channel characteristics. 

30 To the purpose, it is introduced a forgetting factor lower than one unit, as indicated on 
the other hand also in the known art, which appropriately weighs the elements of some 
matrices storing the historical data used in the estimate, before any updating iteration 
of the same. 

The equalization degree obtained with the method described is generally 
35 sufficient to enable a correct demodulation in almost all the real transmission 
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situations. However, in presence of high-intensity non-stationary interferers 
[C 1 1 <-\QdB), wrong decisions taken during the starting phase of the adaptation, 

immediately after the arrival of the noise, are so numerous to overcome the threshold 
enabling the estimate of the best equalization parameters to converge towards BER 
5 acceptable values in the demodulated signal, in a reasonably short time versus the 
semiburst length. The high increase of the error probability is unavoidable due to the 
intrinsic delay of short term decisions, which slackens the instant of interferer view, and 
that of the implementation of the countermeasures to neutralize the same. 

Characteristic part of the claim 1 discloses the way to face high intensity 

10 non-stationary interferers. Supporting matter in the text describes an embodiment 
including a phase of continuous monitoring of the power received on the antennas, and 
of non-coherent detection of possible peaks higher than a given threshold. This phase 
is profitably performed in parallel to the normal operation. The detection of a peak at an 
instant n is used to automatically select the operation according to the embodiment. 

15 After that, for the whole AT duration (typically twenty samples) of a convergence 
transient of the matrices of historical data used in the estimate, only the coefficients of 
the spatial filter are updated and not the estimated samples of the channel impulse 
response, which maintain the values they had at instant n . This strategy is aimed at 
anticipating at instant n the future knowledge on the interferer to cancel, present in the 

20 samples of the stored burst. In this way the interferer is "seen" in advance in the 
calculations and the adaptive error minimization algorithm can start updating the 
beamformer coefficients to try to eliminate the noise. Of course such an advance is not 
possible for the samples of the impulse response, whose adaptive estimate 
mandatorily requires short term decided data. 

25 The estimate as in the characteristic embodiment results therefore Vitiated", 

because to obtain the advantages it aims at, it introduces a synchronization loss 
between the samples of the signals coming from the antennas and the relevant 
decided data. The choice to maintain the estimated samples of the channel impulse 
response constant, reveals therefore as a caution measure destined to avoid error 

30 propagation, otherwise due to the sudden updating of the sole spatial component of the 
structured matrices of statistical nature. 

The original operation is automatically re-seiected when the transient of AT 
duration elapses. 

The advantage obtained by the characteristic embodiment is to render the 
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convergence of the equalization parameters estimate more effective, also in presence 
of sudden and strong interferers enabling out of the midamble, reaching at the end of 
the brief transient a more robust cancellation of the echoes of the interferes 

To make the point on the sequence of equalization phases implemented by the 
5 subject method of the invention for the demodulation of the whole burst, the following 
sequence is obtained: equalization starts in presence of the training sequence, once 
the midamble is exhausted it proceeds using low delay estimated data, in the 
meanwhile the power on the antennas is continuously checked, if before the end of the 
burst the power is over threshold enters the characteristic embodiment of the invention 

10 for rejecting the strong asynchronous interferent, the embodiment quickly exhausts its 
operation within the time AT, at the end of which the previous operation is 
automatically reselected, all within the time interval of the 0,577 ms burst. 

The method of the invention performs the total filtering in two subsequent 
phases, a first only spatial phase, and a second only temporal, but to optimize at the 

15 same time the coefficients of both the filters (as for the coefficients of the time matched 
filter, the same are indirectly optimized through optimization of the samples of the 
impulse response). Thus a reciprocal iteration determines between the two types of 
filtering tending to the optimum, due to the effect of the constraint imposed by the 
MMSE criterion (Minimun-MSE) simultaneously on the two filters. The result of such an 

20 interaction is to aid the crossing through the cascade of the two filters of the maximum 
possible power coming from the useful signal, therefore including its echoes and to 
reject the major possible number of co-channel interferers, to which the relevant 
echoes must be related. What is actually made, is an automatic ^synchronization of 
echoes, valid both for the useful and for the interferers, which operates in synergy with 

25 the discrimination of the direction of the different signals on the azimuth plane. 

The term "joint estimate" or "joint optimization" in the present application 
assumes an effective meaning that it shall be appreciate when disclosing the 
mathematical algorithm of optimization, either in closed or in iterative form. Concerning 
the closed form, the method obtains a system of algebraic equation whose solution 

30 furnishes a vector of the optimum weights and a vector of the optimum samples of the 
channel impulse response. The method consists of the systematic application of other 
methods (Lagrange multipliers, calculation of the Schur complement of a matrix R, 
estimate of an eigenvector associated to the minimum module eigenvalue of the Schur 
complement) each known per s6, but the fact shouldn't prevent the originality of the 

35 Applicant's combination in the particular applicative context. The disclosed solution 
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really minimizes the objective function, maximizing in the meanwhile the Signal to 
Noise+lnterference Ratio (SINR) by mixing in the matrix R, spatial, temporal, and 
spatial/temporal covariance matrices. Consequently the concept of "joint" extends up to 
the finding of solutions. In the disclosed approach the space weights and time samples 
5 are searched simultaneously by setting up a system of equations, each equation is the 
gradient of the objective function, this system gives the "joint optimization". The solution 
of the system is far from being linear as the space and time solutions are inter-related. 
However, the system of equations belongs to the class of separable system that has a 
closed form solution. The inter-relationship between the two spatial and temporal terms 

10 gives the optimum solution as the eigenvector associated to the minimum module 
eigenvalue of the Schur complement. Updating is carried out by first updating the 
spatial covariance matrix R w and then by inverting this matrix (it is a Direct Matrix 
Inversion (DMI) method where the updating follows the matrix R w estimation). The 
inter-relationships generate synergies that improve equalization resolution and 

15 speed-up convergence, this fact is widely proved in the various scenarios depicted in 
the relevant Figures. 

The alternative approach implemented in the Applicant's invention, namely the 
updating method in iterative form, follows the displacement theory of Cholesky 
algorithm, adapted for space-time processing according to the proposed objective 

20 function. Matrix R, being positive defined, assures the existence of a unique Cholesky 
factor L , consequently the inter-relationship between the spatial and temporal terms 
of the matrix R gives the optimum solution even in the Cholesky context. 

In the RADEMACHER's method the optimization involves three disjoint 
gradients, one independently from the others (and not simultaneously in a system as in 

25 the method of the present invention) by using the same error (probably this originates 
the term "joint"); inter-relationships are neglected and consequent synergies that 
improve equalization resolution and speed-up convergence are lost. 

The equalizer based on the subject method does not simply repropose a device 
described in a functional way retracing the steps of the method, but it exploits the 

30 capabilities of systolic arrays implemented through digital integrated circuits of the 
ASIC type (Application Specific Integrated Circuit) purposely configured. 

Consequently, additional object of the present invention is an equalizer operating 
according to the equalization method of the invention, including systolic arrays for the 
combined estimate of the optimum equalization parameters, as better described in the 
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relevant dependent claims. 

Just for the seek of completeness in rebut the prior art, an ultimate document is 
analysed, namely the European patent application EP-A-0793359, inventors 
MORELAND et all, titled: "Method and apparatus for block decision feedback 
5 equalisation using a Cholesky algorithm". An independent method claim 1 recites 
textually: U A block decision feedback equalizer comprising: 

- an input for a received data signal, the received data signal having a known probe 
sequence and a block of data symbols; 

- a channel impulse response estimator coupled to the input and having a channel 
10 history generator responsive to the known probe sequence; 

- a block detector coupled to the input and the channel impulse estimator for 
detecting each data symbol of the block; and 

- a probe data sequence generator coupled to the channel impulse response 
estimator and the block detector. 

15 MORELAND's invention discloses an improved block decision feedback 

equalizer using the Cholesky algorithm in different embodiments listed at page 3 of the 
cited document. Among those embodiments the unique in such a way comparable with 
the equalizer of the present invention is the one in which the block detector includes a 
diversity combiner and a time varying channel estimator. The equalizer of the citation is 

20 designed exclusively for a DFE structure that is know be quite different from an 
equalizer based on a Viterbi estimator. The main difference between a DFE equalizer 
and an equalizer based on a Viterbi estimator (named for sake of brevity Viterbi 
equalizer) is that the first makes a hard decision symbol-per-symbol, while the second 
a soft decision based on the choice of a sequence reproducing a survivor path in a 

25 trellis, after the corresponding latency. Outside the midamble a hypothetical model of 
time-variant channel must be introduced, due to the impossibility to have a trellis and a 
store of the relevant paths for short-time decisions. Concerning the use of cited 
equalizer in presence of time-variant channels due to echo Doppler and/or multipath 
propagation, the best proposed solution disclosed in the document is a block detector 

30 including a space diversity combiner and a time-variant channel estimator. As 
universally known a space diversity combiner, in the most practical cases, includes two 
opportunely spaced antennas respectively connected to a main and a diversity 
receiver. The main and diversity received signals are summed up after the reciprocal 
phase offset and, in case, individual gain adjustment in order to increase the quality of 
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the combined signal. Even combining more than two received signals, the diversity 
combiner acts in a totally different way from a beamformer. In fact the reciprocal phase 
offset introduced in diversity combining is generally random, on the contrary of a 
combination subject to a beamforming algorithm in which, starting from a condition of 
5 fixed phase offset between all adjacent antennas, each subsequent reciprocal phase 
offset is algorithmically controlled in order to change the initial radiation pattern in a 
way to advantage utiles and reject interferents. The absence of a calculated phase 
offset makes a diversity combiner unsuitable for beamforming and consequently for the 
employment in space-time equalizers. In conclusion, a multiplicity of diversity antennas 

10 is not a directive array, this fact has a counterpart in a different structure of the spatial 
covariance matrix and the data covariance matrix. In addition, the absence of spatial 
and temporal interrelation in the solution of the minimization algorithm, deprives the 
DFE equalizer of the cross-covariance matrix between data and inputs to the antennas. 
For these grounds the Cholesky structure of the matrices depicted in matrix equation 

15 45 of the citation is different from that used by the Applicant to introduce a global 
Cholesky factor, as cleared up later. 

The outlined differences between the equalizer of the citation and the adaptive 
equalizer of the present invention, both in the architecture and operation, are so 
relevant that this citation, although disclosing a possible use of the Cholesky algorithm 

20 in an equalizer, discloses an use unfeasible in a space-time equalizer like the 
Applicant's one. 

BRIEF DESCRIPTION OF DRAWINGS 
The many advantages and scopes of the present invention will be appreciated 
from the following detailed description of a preferred embodiment of the same with 
25 reference to the attached drawings, in which: 

- Fig. 1 shows an array of antennas (array) used in the present invention; 

- Figures 2 and 3 show the block diagrams of a same number of beamformers made 
according to the prior art; 

- Figures 4 and 5 are the radiation diagrams on the azimuth plane, obtained applying 
30 the beamformer of fig. 3 to an array similar to that of fig. 1 , concerning two different 

transmission scenarios in mobile radio environment; 

- Fig. 6 is a general block diagram of a receiver particularly emphasizing an equalizer 
operating according to the equalization method of the present invention; 

- Figures 7 and 8 are two radiation diagrams on the azimuth plane, obtained applying 
35 to the array of fig. 1 a beamformer included in the equalizer of fig. 6, relating to two 
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transmission scenarios identical to those that originated the radiation diagrams of 
figures 4 and 5; 

- Fig. 9 shows two time successions of images placed side by side, ordered from top 
to bottom for growing values of the discrete instant n, out of which, the left 

5 succession indicates the radiation diagrams on the azimuth plane obtained applying 
the beamformer of fig. 6 to the array of fig. 1 . and the right one shows the estimated 
channel impulse response; 

- figures 10 and 11 show the BER curves at the output of the receiver of fig. 6, 
according to the possible length of the training sequence of the transmission burst, 

10 for two different values of C/7; 

- Fig. 12 shows a general block diagram of a receiver particularly highlighting an 
equalizer operating according to the equalization method of the present invention; 

- Fig. 13 shows a trellis relevant to a STIM block included in equalizers of figures 6 
and 12; 

15 - Fig. 14 shows two trends of BER along the burst, evaluated at the outputs of 
receivers in figure 6 and 12, respectively; 

- Figures 15, 16, and 17 show the radiation diagrams on the azimuth plane relevant to 
the curve shown with solid line in fig. 14; 

- Fig. 18 shows three trends of BER along the burst, evaluated at the output of 
20 receiver of fig. 1 2, for three different values of the C/l ratio of the input signal; 

- Fig. 19 shows to comparison purposes two trends of BER along the burst, evaluated 
at the output of receiver of fig. 1 2; 

- Fig. 20 is a block diagram of an estimate block FUNZ included an equalizers of 
figures 6 and 12 of the invention; 

25 - Fig. 21 is a block representation of a first systolic array included in the FUNZ block 
of fig. 20; 

- fig. 22 is a block representation of a second systolic array included in the FUNZ 
block of fig. 20; 

- fig. 23 shows a class of additional blocks included in the second systolic array of fig. 
30 22; and 

- fig. 24 shows the signals involving the generic blocks of the systolic array of figures 
21 and 22. 

DETAILED DESCRIPTION 
Making reference to fig. 6, we notice the general block diagram of a receiver 
35 adopting an equalizer operating according to the method of the invention. The receiver 



DB 755 PCT 




is connected to an array ARY of M antennas giving signals, represented as a whole 
with an algebraic vector y(w). The vector y(n) reaches respectively the inputs of a 
block BMF and the first inputs of a block FUNZ. The BMF block is a spatial filter, or 
beamformer, reached by a vector of the optimized weights w op ,(n), coming from first 
5 outputs of the FUNZ block. At the output of the BMF block there are the real and 
imaginary (l,Q) components of a signal y out {n) which is filtered by a FADAT filter 
matched to the impulse response of the transmission channel. The. components of the 
base band signal filtered by FADAT reach the first inputs of an estimator STIM, at the 

output of which an estimated sequence {x a } is present. The second input of the FUNZ 

10 block is reached by a training sequence {x n } stored in the receiver, previously 
synchronized to a copy of the same present in the midamble of the transmission burst, 
availing to the purpose, known methods based on the oversampling of the signal 
received and the correlation of the sequences generated by the oversampling with the 
sequence stored in the receiver. A vector of L optimized samples h opr (n) of the 

15 channel impulse response comes out from second outputs of the FUNZ block, directed 
towards additional inputs of the FADAT filter and towards second inputs of the STIM 
block. 

The FUNZ block includes in its turn a PESI block that calculates a vector of M 
weights w(n) of the beamformer BMF, and a CAMP block that calculates a vector 

20 h(w) of L elements, corresponding to a same number of significant samples of the 
impulse response of the transmission channel. It includes also a two-input adder 
£ reached by two digital signals obtained from appropriate processing of PESI and 
CAMP blocks, from which the signal coming from CAMP is taken with sign minus. At 
the output of the adder an error signal e n is present that returns to the above 

25 mentioned blocks. The first and second inputs of the FUNZ block correspond therefore 
to the inputs of PES! and CAMP filters, respectively, the outputs of the same filters are 
also the first and second outputs of the above mentioned block. 

The representation of the FUNZ block highlights only the main aspects of the 
calculation of w 0/7 ,(n) and h op( {n), a more detailed description of the internal 

30 architecture of this block shall be made describing fig. 20. 

For ease of discussion, the diagram of fig. 6 does not include some blocks 
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referred to known operations made by the receiver and necessary to the operation of 
the equalizer, such as for instance those relevant to the carrier synchronization, the 
conversion/demodulation in base band of the M bursts received from a same number 
of antenna, the sampling and analog/digital conversion of the same, the storing of the 
5 M bursts of digitized symbols, the acknowledgement of the training sequence and the 
synchronization of stored bursts, etc. Furthermore, all the elements of vectors are 
intended referred to complex envelopes in base band. In particular, the y(n) signal 
coming from the antenna is that defined in the (21). 

In the operation, the spatial filter BMF is a MTAPS transversal filter obtaining the 
10 signal y^in) at its output making the weighed sum of the y(n) signals coming from 
the antennas and present in the bursts stored: 

^(») = w r y(«) (13) 
The adapted filter FAD AT and the estimator STIM require, for their correct 
operation, to know the vector h opt {n) t after that their realization is known to the 

15 technician of the sector. In particular, the FADAT filter is a LTAPS transversal filter, in 
general with frequency response of root raised cosine type, where L is the channel 
length. On the contrary, the STIM block is a processor performing an MLSE estimate of 
the transmitted sequence on the basis of the well known Viterbi algorithm, applied to a 
trellis with branch metrics whose expression is the known one of the family of receivers 

20 operating with pilot symbols (training sequence) and differential coding. 

In the equalization method governing the operation of the equalizer of fig. 6, the 
estimate of vectors w^, (n) and h opt (n) is made jointly, considering the following cost 
function: 

r 2i M 1-1 



l=l i=0 



(27) 



25 from which it results that the mean square error e 2 , differently from the Wiener-Hopf 
optimum combiner, results a function both of coefficients w, of the beamforrner and of 
the samples of the channel impulse response h t . In this case the time dimension is 
directly introduced in the error functional to minimize. 

It is now proposed a matrix notation, to the purpose to render the processing of 

30 the problem of the minimum of the (27) more compact. Denote with Y the matrix 
corresponding to N time samples at the M antennas: 
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Y = 



y r (n-l) 



(28) 



_y r (n-Af+ l)_ 

obtained from vector y(n) = [>>, (/?), y 2 («), y M (")] T whose M elements, are the 

samples of signals received on the antennas; therefore: 



Y = 



y n -\ 



y n 



y n *M-\ 
y n +M-i 



(29) 



Concerning the number we shall see in the description of figures 10 and 11 
that a value included between 20 and 30 ( N = 26 for the GSM midamble) is sufficient 
to furnish a good statistical basis for the estimate of the best equalization parameters. 

Denote with X the matrix of data transmitted (limited just for now to the training 
sequence), each one having duration L: 

x r ("-l) 



X = 



(30) 



_x T {n-N + l)_ 

obtained from the vector x(n) = [x n (n) 9 Jc lf+l (nX....,x -l+I . l (/i)] r whose L elements are 

data transmitted coinciding with the training sequence {x n } , and duly stored, assuming 
a channel length equal to L\ so that: 



X = 



x 



n-l 



v n-Af+2 



n-N+L . 



(31) 



15 and with: 



w 



2 »■ 



h = [/*,, Z^, 

the vectors of the coefficients of the BMF spatial filter and of the channel impulse 
response are respectively indicated. 



•A] r 



(32) 
(33) 
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The problem proposed in the equation (27) can therefore be rewritten in the 
matrix form, as follows: 

( w om> h op t ) = argmine 2 (w,h) = arg min||Yw - Xh|| 2 (34) 

\ y H * w.h .h 

l|h|| 2 = 1 (35) 
5 where in the (35) is introduced the constraint on the power of the impulse response to 
avoid the banal solution (w^ =0, h opt =o) . The constraint (35) is a consequence of 

the fact that the bits of the sequence transmitted are assumed as time uncorrelate, 
being valid other hypothesis the (35) should be modified in a consistent way. 

The problem of the minimum constrained can be faced applying the Lagrange 
10 method. To the purpose, denoted with y a real parameter (Lagrange multiplier), let's 
consider the following functional: 

J(w, h, A) = || Yw - Xh|| 2 + r (llhll 2 - l) 06) 

The research of the minimum constrained in the (34) reduces then to the 
research of a "solution" (w, h) of the problem of free minimum of the functional 
15 j(yv,h,y) . The above mentioned solution can be searched considering the system: 

J w (w,h,y) = 2Y*(Yw-Xh) = 0 (37) 
y h (w,h,x) = -2X*(Yw-Xh) + 2>* = 0 (38) 

^(w,h,r) = h'h-1 = 0 (39) 

The equation (37) can be construed as the introduction of the error orthogonality 
20 principle, in fact it is possible to express the equation (37) as: 

4 e -y'(")]=° (40) 

Resuming, for convenience, some expressions partly already introduced, denote 
with the spatial covariance matrix expressed by: 

R„ = E[y(t)y'(t)] (41) 
25 and likewise with R a and R v , respectively, the covariance matrices on the data and 
of cross-covariance between data and inputs to the antennas expressed as: 



aj 

3J_ 
do. 

aj_ 
ax 




= E[x{t)x{t)] (42) 

R„ = £[x(r)y'(/)] (43) 
We shall use below the sampling co-variance matrices defined as follows: 

^ n-l 

R«=^2>(«K(") <"> 

R v =^f>(«)y"(») (45) 



We can therefore define a matrix R structured as follows: 

R ... R' 



R 



yy « XV 

R_ R^ 



(46) 



which shall enable a more simple representation of the problem. Equations (37) and 
10 (38) can be re-written respectively as: 

R^w = R\,h (47) 
R^w-R^h = yh (48) 

From the equation (47) it is possible to obtain w as function of h and replace it 
in the (27) to obtain the mean square error. While replacing the w value in the (48) the 
15 following equality can be demonstrated: 

e 2 (w,h) = h-R^h = - r ( 49 ) 

where: 

(so) 

is called Schur complement of R^ in R. The estimate of vector h(n) that minimizes 
20 the mean square error, corresponds therefore to the estimate of the eigenvector 
h (n) of matrix R 1 associated to the minimum module eigenvalue y . 

opt * * s 

Therefore it is possible to estimate also the w opt (n) coefficients of the 
beamformer replacing vector h opt {n) in (47): 

25 Summarizing, the choice of the optimum weights minimizing the cost function, 
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can be synthesised in the following main points: 

1 . calculation of R , matrices and of the Schur complement Rf ; 

2. estimate of the eigenvector h opt (n) of matrix Rf associated to the minimum 
module eigenvalue / ; 

5 3. estimate of w 0/7 ,(h) coefficients of the beamformer. 

Assuming to implement the equalization method summarized under points 1., 2., 
3 M with a program controlling a finite precision processor, the calculation of the Schur 
complement R^ can result a very critical operation, mainly if the word length of the 

processor is very limited. 

10 The above mentioned criticality is considerably reduced thanks to the use of the 

methods developed in a theory on partial triangular factorization of the matrices, known 
in literature as "Displacement Theory", or even "Displacement Structure", included in 
the bibliography given below, reference to which shall be made. The key element of 
this theory is an entity called "Choiesky factor" of a matrix, in the case considered 

15 matrix R, whose calculation is made through the simple updating of a previous value, 
intervening with appropriate transformations (Givens rotations) directly on the new 
data. As we shall see, the updating of the "Choiesky factor" gives a matrix, structurally 
different from the R but having the same information content, whose elements highly 
facilitate the w op , and h opt calculation. It is thus excluded the explicit calculation of 

2p sampling covariance matrices and the recourse to critical operations, such as that of 
the direct matrix inversion, since the above mentioned rotations do not involve any 
calculation of the determinant. This procedure involves a better conditioning of the 
calculation problem, and a consequent reduction of the minimum word length of the 
processor. 

25 Such an approach gives the additional advantage to extend the calculation to the 

whole burst; the solution described is in fact intrinsically recursive (iterative) . On the 
other hand, a non recursive method, that is involving the explicit calculation of the 
sampling covariance matrices, is characterized by a high complexity in the calculation 
of the Schur complement Rj through the (50), which makes the extension of such a 

30 method to the whole burst difficult. 

It must be underlined that the application examples of the "Displacement Theory" 
known up to now, concern only the optimization of the performances of time-only 
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equalizers, or time equalizers used in space-diversity combiners as in the cited 
MORELAND's invention. A new element has been that to connect the mathematical 
formulation of the solution of the problem of minimum of the error functional used in the 
present invention, to the formulation given by such a theory for the solution of the 
5 problem of minimum in MMSE-DFE time-only equalizers. Due to the similarity 
highlighted between the two formulations it has been therefore possible to avail of all 
the knowledge . developed for time equalizers adjusting it to the space and time 
equalizer of the invention. 

This stated, let's consider the matrix R described in the (46) which, being 

10 positive defined, assures the existence of a unique Cholesky factor L such that: 

R = EF (52) 
One of the fundamental methods for the triangular factorization, is the Schur 
reduction process described: 

• in the article by J. Chun and T. Kailath, under the title "Generalized Displacement 
15 Structure for Block-Toeplitz, Toeplitz-Block, and Toeplitz-derived Matrices", 

published in NATO ASI Series, vol. F 70; 

• in the volume "Numerical Linear Algebra Digital Signal Processing and Parallel 
Algorithms", authors G. H. Golub and P Van Dooren, edited by Springer-Vetlag, 
1991. 

20 The process enables, after M factorization steps, the reduction of matrix R as 

follows: 



R = 



L 

U 



[L- U-] 



0 0 

0 R^ 



(53) 



where Rf can be then considered as the residue of the partial displacement in 

triangular form. This approach enables to calculate the matrix R^ without having to 

25 invert the spatial covariance matrix R^ . Starting from the Schur complement, obtained 

from the factorization process, it is possible to estimate h opt according to classical 

estimation methods of the eigenvector associated to the minimum module eigenvalue, 
for instance, the method of inverse powers. 

The partial factorization highlights also the possibility to have matrices R^ and 

30 R^ factorized in their triangular form: 
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=LL' (54) 

= UL' (55) 
enabling the reduction of the equation (51) in the form: 

L"w op , = U'h opJ (56) 

5 The lower triangular form of matrix L, as we shall see below, enables the 

solution of the equation (56) with a simple process called "back-substitution" process. 

Let's now analyse the different sequential steps of the equalization method 
employing the Cholesky factor in the calculations, giving for each of then the relevant 
complexity: 

10 • step 1: application of the known Schur algorithm for the quick reduction of the 
matrix R for the Cholesky factorization; 
• step 2: application of an algorithm known in the literature as "Method of inverse 

powers" for the estimate of the eigenvector {h opt ) of R^ associated to the minimum 

module eigenvalue; 

15 • step 3: solution of the equation (56) for the estimate of optimum coefficients of the 
beamformer (w opt ), according to a method known in literature as "back-substitution" 

method. 
Stepl 

From the Toeplitz theory of matrices it is known that a generic matrix R(/?) has a 
20 displacement structure Toeplitz-type if the displacement matrix VR(«) defined as: 

VR(w) = R(w) - F(/i)R(« - 1)F*(«) (57) 
has low rank (rank(VR)"rank(R) ). The matrix R expressed in the equation (46) falls 
within this category and, in accordance with the above mentioned theory, it can be 
calculated in an iterative way, arbitrarily assuming R(0) = 0 t and extending to the 
25 whole midamble the updating defined as follows: 

R(n) = R(«-l) + g„ g ; (58) 
where the updating vector g n , denoted as: 

'y(n)' 



g„ = 



x(») 



(59) 



is called generator of the difference matrix, since imposing F(«) = 1 we can write 
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VR(n) = g„g n 

Likewise, it is possible to update the Cholesky factor of matrix R(/i), thanks to 
the special time-varying nature of the same. Therefore, denoted with L(n) the 
Cholesky factor of matrix R(n), such that R(n) = L(«)IT(«) , it is possible to verify 
5 that: 

L(n)i7W = [L( W -l) gJ^V 0 

. L 8- - 

This suggests an effective method for the updating of the Cholesky factor 
through the use of unitary transformations, so that, starting from the factor calculated at 
instant n-1 and generator at instant n, it is possible to determine the new factor at 
10 instant n without having to calculate the matrix R(n) . 

In fact, it is possible, arbitrarily assumed L(0) = 0, to express the updated factor 

as: 

[LW 0] = g fl ]Qto (61) 

where Q{n) is a unitary rotation, in general made as a sequence of elementary 
15 transformations, Q 0 ,Q P ' which, applied to [l(w-1) g n ] produces the nulling of 
the column of the generator g„ introducing a zero for each elementary transformation. 

An example of recursive updating, let l Q (n) and / 0 (w-l) the first elements in 
position (0,0) of matrices L(/i) and L(n - l) , respectively and g n 0 the first element of 
the generator vector. It derives that the first line of matrix L(«) can be written as: 

[/„(*) 0 0] = [/ 0 («-l) 0 ^ 0 ]Q 0 (62) 
Applying the same transformation to the whole matrix we obtain: 

0 0 



[L(»-l) g„]Q 0 = 



loW L.0.-1) g„,. 



(63) 



where L ,(«-!) and g rt l are the matrices obtained after elimination of the first 
element of the generator column. 
25 If we recursively proceed in the updating of the whole factor L(«) , we can see it 

factorized in the following triangular form: 
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L(«) = 



(64) 



'L(w) 0 
U(«) L,W_ 

where L and U are the matrices defined in the (53), and L, is the Cholesky factor of 
R^, such for which: L,L* =R S J ". As it can be noticed from the (64), matrices L and 
L i are triangular while the U is necessarily rectangular, or possibly squared. 
5 The factorization of the (64) results being essential since it enables the 

calculation of the vector of weights w op ,(n) without making the inversion of the matrix 

(see equation 50) and, as we shall see in short, without inverting the matrix Rf , 
being this operation requested by the algorithm of inverse powers for the calculation of 
the eigenvector (h opt ) associated to the minimum eigenvalue. 
10 The recursive factorization for the updating of the whole Cholesky factor requires 

a number of complex multiplication equal to 2(M + L)(M + L - 1 ) and a number of sums 
equal to (M+ L)(M + L - 1). 
Step 2 

The method of inverse powers for the calculation of the eigenvector associated 
15 to the minimum eigenvalue requires the iterative solution of the equation: 

K = (R^W)"'h w (»-D (65) 

■vM - & « 66 > 

where the normalization of the (66) is introduced as the result of the constraint on 
power of the impulse response of the estimated channel: |h| 2 = 1. Due to the above 
20 mentioned constraint, it is convenient to select an initial vector h op( (0) with at least one 
non null element. 

The Cholesky factorization enables the calculation of the above mentioned 
eigenvector without having recourse to the inversion of the matrix, since the (65) can 
be re-written as: 

25 L f (n)L»h. ■ h^(»-l) (67) 

The equation (67) can be factorized in two elementary phases, introducing an 
intermediate vector u„ defined as: 

L»h„ = u„ (68) 
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in order to be able to express u n as solution of: 

M»K = M"- 1 ) ( 69 ) 

We can solve the equation (69) regarding u„ and then, having found u n , we can 
solve the equation (68) regarding h n . Both the equations are simple to solve since the 
concerned matrices are triangular matrices. Should we have to solve equations with 
generic square matrices L x L, we should expect a number of arithmetical operation in 
the range of L 3 . On the contrary, being the matrices triangular it is possible to reach the 
solution in only L 2 operations. To illustrate this approach, known in literature and called 
back-substitution, we develop equations (69) and (68): 

l 2i u x +1^ = h opt2 (n-~l) 



(70) 



yIl*«,L +/ I-lA.I-l + "" B+ A*A,l " U l 



(71) 



25 



Starting from the first term of the (70) we can calculate u x : 

»i = ] 

Ml 

and replace it in the other L - 1 equations. This demonstrates that we can reduce the 
problem to a triangular system of dimensions (L - 1) x (L - 1) at the cost of L - 1 
multiplication-and-sum operations and one division. With this method we can solve the 
system at the cost of L divisions and (L - 1) + (L - 2) + ... + 1 multiplication-and-sum 
operation, for a total of !4L(L - 1) multiplication-and-sum operations and L divisions. 
After having solved the equation (69) regarding u n we can solve the equation (68). 
The number of operations required is fully identical to that required for the solution of 
the first triangular system of equations. 
Step 3 

Also the calculation of the optimum coefficients of the beamformer can be made 
for a back-substitution in (56), directly in the space of the signals (without explicit 
calculation of the covariance matrix), calculating: 
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L>)w op >) = V m (*)ti m (n) (56) 

In this case the previous step 2 must b extended introducing also the calculation 
of the second term of the (56), for a total number of ML + 14M2 multiplication-and-sum 
operations and M divisions. 

The following table summarizes data relevant to the calculation complexity in the 
three steps described, giving also the calculation power of an hypothetical processor 
(expressed in MIPS) requested for the real time execution. It is considered an updating 
of coefficients, extended to the sole midamble, for an array of antennas of M = 8 
receivers and for a matched filter having length L = 5. 



Step no. 


No. of complex multiplications 


No. of complex sums 


MIPS 


1 


2(M + L)(M + L- 1) 


{M+L)(M + L- 1) 


8,97 


2 


L(L-1) + 2L 


UL-V 


0,86 


3 


ML + ViM 1 


ML + VzM 2 


0,098 



The imporovement introduced by a joint estimate can be demonstrated 
considering the previous scenario proposed in fig. 4, where the time dispersion 
introduced by the channel is represented by an echo of the signal delayed by 5 /xs. 

Fig. 7 shows the radiation diagram on the azimuth plane determined by the 
beamformer BMF of fig. 6, where it is clear the different behaviour versus the situation 
of fig. 4 concerning the delayed echo. The radiation diagram of fig. 7 shows in fact a 
power level in the echo direction, of the same extent of the level in the useful direction. 
This enables, to the price of a controlled intersymbols interference, to maintain high the 
power of the useful signal present in the estimate process of the transmitted sequence 

The robustness of the new functional proposed can be demonstrated in fig. 8, 
just considering the behaviour of the beamformer BMF in the same critical scenario of 
fig. 5, where a replica of the co-channel interferer was present in the same direction of 
the useful signal. 

As it can be noticed from the radiation diagram of fig. 8, the power is now mainly 
transferred in the useful echo direction, so that the diagram shows a zero also in the 
useful direction, in addition to that of the interferer. This behaviour can be explained 
since the introduction of the time dimension in the functional enables the combination 
of echoes, even when these should be time delayed. We can notice the capability of 
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the mixed functional to select in space, among the different delayed echoes that are 
coherent with the training sequence, the one characterized by the greater ratio 
between signal and interferer. 

This behaviour enables to manage also mobile users with superimposed nominal 
5 directions, whenever they have an opportune spatial dispersion of echoes. It can be 
experimentally demonstrated how, in RCS environment, this enables a considerable 
reduction of the intercell handover probability. 

The joint estimate of the optimum equalization parameters, however requires a 
higher number of parameters than a time only estimate. In particular, in the case of 
10 FDMA-TDMA systems with comparatively short training sequences (just as in the 
GSM), two potential problems arise: 

• a growth in calculation complexity in the joint equalization method, compared to a 
disjoint equalization of the prior art; 

• a higher variance in the estimate of the optimum equalization parameters, and 
15 therefore a slower convergence speed of the same. 

The length of the training sequence represents a critical parameter in the 
equalizers used, or that can be used, in the above mentioned systems, since the 
algorithms for the joint estimate of the optimum equalization parameters must identify 
the channel impulse response, filtered by the beamformer BMF, availing of a limited 
20 number of samples, that is the 26 of the midamble of the GSM burst. 

Fig. 9 shows, from top to bottom, a sequence of radiation diagrams on the 
azimuth plane and of representations of the channel impulse response, "photographed" 
on completion of n iterations of their joint estimate, made through the algorithms 
adopted at steps 1, 2 and 3. The scenario shown relates to a useful signal and to a co- 
25 channel interferer (denoted with u o" and "x n respectively on the outer arc), both 
accompanied by their own echo delayed by 10 (denoted with the same symbols 
placed just under the outer arc). 

Making reference to fig. 9 it is possible to notice how the total convergence is 
dominated by the convergence speed of the impulse response h opt (n). At instant n = 

30 6, when the statistical basis is not sufficiently extended to make a correct estimate of 
the impulse response, the beamformer BMF introduces some zeroes also coinciding 
with the echoes of the useful signal. Also for n = 16, that is when the estimate of 
h opt (n) leads to recognize the sole main impulse, the beamformer BMF does not 
behave yet in an optimum way, since it introduces a zero in a position coinciding with 




the delayed echo of the signal. However, we can already notice that for n = 26, as well 
as for n = 30, the channel identification enables the identification of the delayed echo 
as useful signal, and therefore the radiation diagram shows the zeroes only coinciding 
with the interferer and the relevant echo. 

A more detailed analysis of the convergence speed can be made evaluating the 
BER at the output of the receiver of fig. 6 as the midamble length varies, when it 
operates in critical conditions, for instance with C/l = - 10 dB. 

Figures 10 and 11 show the results of this analysis referred to a scenario that 
includes a co-channel interferer, in addition to the useful, for three different values of 
the signal-noise ratio SNR at the input of the receiver: 10, 20, and 30 dB. More in 
particular, the curves have the length of the training sequence on the x-axis, expressed 
as number of samples, and the BER on the y-axis. Fig. 10 is referred to a situation with 
C/l = - 10 dB, while fig. 1 1 to a C/l - 0 dB. 

Making reference to figures 10 and 11, the receiver performance can be 
considered stabilized in both the cases for values of midamble length between 20 and 
26 samples. In fact, the almost total slope variation of the curves is in this area, and 
this indicates that the convergence of the optimum equalization parameters inside the 
training sequence is reached. 

Making reference to fig. 12, we notice that the general block diagram of a 
receiver whose equalizer is the one shown in fig. 6 with the addition of a portion that 
considers the adaptive operation out of midamble. 

Comparing figures 6 and 12, it can be noticed that the two representations of the 
equalizer differ due to the fact that in fig. 12: 

- and additional RIT block is highlighted, introducing a delay D on the path of the bus 
of the input vector y(n) towards the FUNZ block; 

- the FUNZ block is reached by a vector x f (n) instead of the training sequence [x n } ; 

there is a new DECTEMP block that stored the vector x'(w), coming out from the 
STIM block. 

As we shall see describing the relevant operation, the schematization of fig. 12 is 
suitable to the operation in presence of channels with time varying characteristics 
and/or in presence of non-stationary interferes, that is enabling out of the training 
sequence, which therefore cannot be identified and cancelled considering the sole 
midamble in the estimate of the optimum parameters of the equalizer. Just for this 
reason the two schematizations of figures 6 and 12 are complementary between them, 



DB 755 PCT 



41 



meaning that during the processing of the midamble the first one is active, and at the 
end, and for the whole remaining duration of the burst the second one is active, in 
practice both the schematizations can be referred to two different operation modes of a 
unique structure of equalizer. 
5 As it is known the Doppler effect is one of the causes that in a mobile radio 

context more contributes to render the characteristics of channels time varying. It 
consists in a small frequency variation of the signal following the movement of the 
mobile, which involves a slow time variation of the amplitude of echoes received. The 
channel thus results time-varying, and even if the variations take place very slowly, can 
10 be useful to try to follow them updating time by time the optimum values of parameters 
h opt and w opt estimated during the midamble. In this way we can improve the 

demodulation quality, since the estimator STIM is made adaptive, since it recalculates 
at each step the branch metrics on the basis of the estimate of h opt , which in its turn is 

affected by the w opt one. 

15 From the operational point of view, it may be convenient to update (according to 

the method described later on) the optimum couple of parameters h opt and w opt 

estimated during the midamble, and simultaneously equalize the signal received, 
starting from the two ends of the midamble and proceeding in parallel in the two 
directions towards the burst ends, previously stored and synchronized. The deriving 

20 higher circuit parallelism enables to avail of a double time for the execution of the 
calculations. This is a significant advantage, considering the high number of single 
calculation elements to perform in a symbol time of 3,69 ^is only, which would be 
doubled using two semibursts in parallel. 

The presence of non-stationary interferers inside the transmission bus can be 

25 explained by the fact that the BS relevant to different cells are not necessarily 
synchronized among them. For instance, in the case of an RCS approach to the 
problem of the increase of the spectral efficiency, we can expect to have co-channel 
interferers in cells very close to that where is the mobile transmitting the useful signal. 
Since it is not possible to synchronise all the BS, it is possible that an interferer disturbs 

30 the transmission only during a part of the burst and therefore it cannot be "identified" 
and cancelled during the transmission of the training sequence. 

The equalizer of fig. 12, thanks to the possibility to update during the whole burst 
the couple of optimum parameters h and w opf estimated during the midamble, is 
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thus enabled to face also the case in which one or more interferes disturb the useful 
signal starting from an instant outside the time slot where the training sequence is 
transmitted. 

In the operation, we can say just from now that, with the details we shall give 
5 below, the algorithm developed by the FUNZ block remains almost unchanged. This 
can also be understood from the great similarity existing between the pictures in figures 
6 and 12, made more evident due to the fact that the RIT block, delaying the vector 
y(n) of a discrete time D can be included in the FUNZ block. 

The adaptive operation requires the knowledge of the sequence of data 
10 transmitted {jc n } . When these are not known, that is when the training sequence is no 

more available, the same must be replaced by data [x n ] decided in the demodulation 

process. It is thus introduced an inevitable delay between the samples of the signals to 
the antennas and the corresponding data of the MLSE estimated sequence, both 
necessary to determine the error. In fact, as it can be noticed in fig. 12, before the data 
15 x n relevant to instant n is available, we must wait both for the delay L introduced by 
the matched filter FADAT, and the delay due to the latency time of the Viterbi algorithm 
necessary for the emission of sufficiently reliable decisions. More in particular, if D is 

the total delay {L + latency of Viterbi), data decided of the sequence {x„ } are available 
only at instant n + D, and therefore the estimate of the best equalization parameters 
20 and v/ op( (n) at instant n can avail of the knowledge of inputs on the antenna 

and of data decided only up to instant n - D t that is the vectors y(n - D) and 
x(n - D) . The RIT block compensates the delay already existing on the other way. 

It is worth to say again that vector x{n) contains now the decided data Jc n , while 
during the combined estimate of the optimum equalization parameters, made during 
25 the midamble, the same vector contained already known data x n and therefore without 
delay. The statement made on vector y(w - D) means that, concerning the calculation 
of h opt (n) and w^(n) , the mean at the n -th instant calculated on the N elements of 
the columns of matrices R^R^ara/R^, concerns elements included between the 

instant n - D and instant n-(D + N). 
30 What in pratice occurs due to the delay D , is that the combined estimate of the 
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optimum equalization parameters at instant n results in delay compared to the channel 
variations, which have effect without any delay on symbols y(n) received. Two needs 
therefore conflict. The first one is to maintain high the latency of the Viterbi algorithm, to 
assure correct decisions avoiding error propagation in the estimate loop. The second 
5 need is to furnish the data decided in the shortest time to the FUNZ block, in order not 
to introduce an excessive delay in the updating of the best equalization parameters, 
such to compromise the possibility to follow quick variations in the channel 
characteristics. 

The approach followed in the subject method is to use in the estimate of the 
10 optimum equalization parameters, a sequence of low delay estimated data, 
^ represented by the L elements of the vector x % {n). The latter, is obtained forcing the 

\D estimator STIM to decide just after few steps ( 4 or 5 approx.), which is the more 

*5 probable sequence, obtained climbing the trellis the maximum likelihood path. Data 

t: finally decided by the sequence {*(«)} are on the contrary available at the end of the 

CO 15 burst, in order the same are as reliable as possible. 

ry The above mentioned technique, which solves the problem to avail in advance of 

Q the data of an almost correctly estimated sequence, is called "Per-Survivor 

12 Processing", or PSP technique, and is described in the paper under the title "Per- 

fln survivor processing: a general approach to MLSE in uncertain environments", authors: 

y 20 R. Raheli, A. Polydoros and C. K. Tzou, published in IEEE Trans. Commun., vol. 43, p. 

354-364, February/March/April 1995. With this technique the Viterbi algorithm can be 
used in presence of one or more unknown, or not well determined parameters, in this 
case the samples of the optimized impulse response, which if on the contrary were 
known could be profitably used in the updating of branch metrics. 
25 Since low delay decisions can result wrong, an error propagation phenomenon 

can be originated due to a vicious circle establishing in the loop formed by the 
feedback of x'(n) toward the block FUNZ. In fact, considering all the elements of 
vector x'(/j) as the result of single low delay decisions at symbol time, as belonging to 
a sequence running in a shift register long L (this is what actually occurs in the 
30 representation of the states according to Viterbi), a wrong decision would finally 
produce its negative effects for a number of symbol times equal to that defined for low 
delay decisions. 

The subject equalization method eliminates also this second drawback by 
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replacing, at any symbol time n, all the L elements of vector x'(w) by a same number 
of temporary data directly copied from the maximum likelihood path in the trellis valid 
at the time of n-th symbol. 

Fig. 13 shows in fact an example where at instant N a wrong temporary decision 
5 is taken in the trellis, which however does not propagate because at the subsequent 
instant N + 1 the decided sequence copied is the correct one. Actually, not always 
things can proceed like this, but very likely this is verified in presence of interferers not 
particularly intense. 

The operation method just described makes the updating of the optimum 
10 parameters out of the midamble more accurate, remaining the approximation 
introduced by low delay decisions. Furthermore, to improve the adaptive behaviour of 
the equalizer of fig. 12, it is necessary to introduce a X coefficient so-called "forgetting 
factor", in order to less weigh the less recent data. By this, recalling the (58), the matrix 
R at instant n can be written as: 
15 RW = ^R(/i-l)4-g n g; (72) 

where X is an arbitrary real constant for which: 0 < A < 1 , and g n is an updating vector 
defined as in the (59), with the difference that its elements are here the low delay 
decisions. 

20 Likewise, recalling the difference matrix (57) and the (58), we can write: 

RW = FWR(,-l)F'W + gflg ; (74) 

being: 

F(«) = VXl (75) 
With this position the updating of the Cholesky factor can be directly obtained 
25 from equation (74). In fact, let L the Cholesky factor of R , we can update the same, 
likewise we have done in the (61), according to the relation: 

[L(n) 0] = [VAL(« - 1) g„]Q(n) (76) 
With all the above clarifications, the estimate of parameters h^(w) and w^n) , 

made adaptive, is obtained applying the recursive algorithm already described for the 
30 non adaptive estimate made by the FUNZ block of fig. 6. 

Making reference to figures 14, 15, 16, and 17, the actual adaptive equalization 



8n = 



(73) 




of the subject method are shown. 

Fig. 14 shows two BER trends along the burst of a same number of demodulated 
signals, for an interferer disturbing up to the 50th symbol. A first trend (dashed line) is 
relevant to the subject equalisation method, assuming to 

freeze the estimate out of the midamble, while the second trend (solid line) relates to 
the same method enabling the performance of the adaptive operation. The input signal 
E C 

shows — - = 10 dB and — = 0 dB ; the trends are averaged on 500 performances. The 

#0 1 

BER indicates the probability that the generic symbol is demodulated in a wrong way, 
for instance for bad equalization. As it can be noticed from the first curve, if the burst is 
demodulated using the equalization parameters estimated with the training sequence, 
the performance greatly deteriorate for the time the interferer is present. 

The second curve shows on the contrary, that the adaptive operation can limit 
the degradation of the performances. In fact, proceeding from centre to the ends, it can 
be noticed how the performances deteriorate coinciding with the arrival (in the 
calculation) of the interferer, to better later on, after equalization. This means that the 
equalizer could adapt to the new transmission situation and vary the radiation diagram 
of the antenna array in order to cancel the interferer. The next figures 15, 16, and 17 
show radiation diagrams on the azimuth plane relevant to the curve with solid line of 
fig. 14. 

Fig. 15 shows the radiation diagram estimated at the end of the midamble, when 
the interferer is not still present. The arrival directions of the echoes of the useful signal 
are represented with solid line with line lengths having length proportional to the power 
of the relevant echo, expressed in dB. 

Fig. 16 shows the evolution of the previous radiation diagram immediately after 
the identification of the interferer, whose echoes are represented by dashed lines 
having length proportional to the power of the relevant echo, expressed in dB. 

Fig. 17 shows the radiation diagram when the joint estimate of the optimum 
equalization parameters has reached the convergence. As it can be noticed, the 
beamformer can cancel the signal coinciding with the arrival directions of the 
interferers. 

The evaluation of figures 14 and 17 shows that applying the equalization method 
scope of the present invention, the joint estimate of the optimum equalization 
parameters can converge towards a couple of vectors h opt {n) and w D/rt (w) that can 
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assure a good quality in demodulation, even in presence of non-stationary interferers. 
What just stated is true until the non-stationary interferer is in the same size range of 
the useful signal. In the contrary instance, a sudden peak of the demodulated signal 
can occur in the BER, that can hardly be taken to acceptable values in a reasonably 
5 short time compared to the residual duration of the burst. 

Fig. 18 shows three trends of BER along the burst of a same number of 
demodulated signals, for an interferer disturbing up to the 50th symbol. The three 

E b 

curves are referred to an input signals with "rr"= 10 dB and three different values of 
C 

— , of 0, - 10, - 20 dB, respectively; the trends are averaged over 500 performances. 

10 The worsening of the convergence as the ratio C/I decreases is evident, in particular, 
we see that the BER does not drop at acceptable values until the ratio C/I remains 
lower than -10 dB. 

The cause for this behaviour can be due to the delay, even small, affecting the 
low delay decisions in the MLSE estimate of data transmitted. 

15 Just due to this delay the equalizer does not immediately detect the arrival of a 

strong interferer, in such a way that, when the first temporary decisions that consider 
the interferer are available, therefore suitable to take measure to cancel the same, it 
will have in the meantime caused a great number of wrong decisions, having an extent 
such to compromise the convergence of the estimate made by the FUNZ block. 

20 The inconvenient is highly attenuated by the original embodiment of the present 

invention, which tries to cancel the delayed vision of the interferer in the adaptive 
equalization, in order to quickly cancel the non-stationary strong interferers. 

The embodiment includes a phase for the continuous monitoring of power 
received on the antennas, and for the detection of possible peaks higher than a set 

25 threshold. The detection of a peak at an instant n is used to automatically select the 
operation of the equalizer according to the methods of the embodiment. Since noises 
to be detected are sufficiently strong ones only, a particularly sophisticated power peak 
detection method is not necessary. For instance, it is possible to use a non-coherent 
detection method, known to those skilled in the art. In the receiver of fig. 12, for ease of 

30 representation, we have not included an appropriate block performing the above 
mentioned functions. 

The detection phase is profitably made in parallel to normal operation out of the 
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midamble. The embodiment quickly exhausts its effects in a transient time having AT 
duration (of about twenty samples) necessary, as we shall see later on t to the setting in 
the computation of some matrices used by the FUNZ block. When AT elapses, it is 
convenient to automatically reselect the original operation that results more effective at 
transient exhausted. 

The strategy implemented by the embodiment to cope with strong intensity 
stationary interferers, is to advance at instant n the future knowledge about the 
interferer to cancel. In this way the interferer is "seen" in advance in the calculations, 
that is before the first relevant decisions are available and the adaptive algorithm of 
error minimization can start updating the beamformer coefficients to try to eliminate the 
noise. The above mentioned strategy is implemented in transient AT , dividing in two 
elementary phases (1 . and 2.) the updating of matrix R(n) . In practice, we proceed as 
follows: 

1. during this first phase, which elapses within the symbol time n , it is anticipated the 
updating of the sole spatial covariance matrix with vectors y(n) obtained from 
the samples of signals on the antennas, present starting from instant n up to instant 
n + AT in the stored burst. The statistical basis of the above mentioned matrix is 
thus extended, that now contains sufficient information of the interferer to cancel. 
The matrix is no more modified up to the end of the transient. Then we proceed 

to phase two. 

2. During this second phase, for n < n < n + AT, the covariance matrices on R^ 
data with vector x'(n) f and the cross-covariance matrices R^ with vectors x'(n) 
and y(«) downstream the RIT block are updated. Concerning vector x'(n) it is 
worthwhile pointing out that starting from n , once the low delay decisions term is 
elapsed and in lack of appropriate measures which shall be shortly specified, this 
vector could express completely wrong decisions due to the performance of phase 
1. 

When the n + AT instant is reached the estimate of the mentioned matrices 
proceeds in the conventional manner. The updating process of matrix R(w), in the time 

interval \n , Ti + AT] , can therefore be expressed as follows: 



R(») = 



(77) 
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expressing the different time interval used for R w calculation with respect to the other 
terms. 

The "anomalous" updating of matrix R(«) generates a setting transient having 
duration AT in the computation of matrices R„ and R^. Said interval just coincides 
with the amplitude of the time shifting in the updating of the spatial covariance matrix 

The joint estimate of the optimum equalization parameters will therefore result 
'Vitiated" by this type of approach, which however proves to be necessary to cope with 
quick variations of the channel characteristics, enabling the beamformer BMF to 
introduce some zeroes coinciding with the echoes of the new interferer, which 
otherwise it will not be able to introduce, unless with an excessive delay . 

As a consequence of the matrix setting time R(«) too many, errors could 
propagate in the estimate loop during the transient [n, n + AT] , therefore a further 
feature of the method embodiment, is to update the vector w 0/>f («) of the weights of 
the beamformer BMF during the transient only, maintaining the vector of the samples of 
the impulse response constant for the whole transient time at a value h opl it had at 
instant n . The choice rewards the coefficients of the beamformer BMF since they 
result less sensitive to decision errors compared to the samples of the impulse 
response. 

For what said above, we can now estimate even during the transient 
[n, n + AT] the optimum equalization parameters. In particular, h op( is constant and 
the (51). recalling the (77), enables to determine the vector w^,, solution of the 
equation: 

Ryy{n + AT)* op ,{n) = R^OOh^ (78) 
in which R^(« + AT) is calculated at instant n and then kept constant for the whole 
AT, while R^(«) is calculated at each Instant w-th of the transient As it can be 
noticed, due to h^, constant, the (67) and (68) need no more to be solved. 

opt 

The previous considerations on the updating methods of matrix R(«) of the (77), 
in the time interval [n, n+AT], are verbatim transferred to the updating of the 




(79) 



Cholesky factor of the same. 

Concerning this last point, at an instant n of the detection of the power peak 

received on the antennas, the following actions are made: 

- within the symbol time n , it is updated the sole triangular matrix L with vectors 
5 y{n) obtained from the samples of signals on the antennas, present in the stored 

burst, starting from instant n up to instant n + AT: 
• for n < n < n + AT , matrices V(n) and L,(«) rectangular and triangular, 
respectively, are updated with said updating vector g n and equation (56) is solved, 
keeping the vector of samples of the impulse response constant at the value 
10 assumed at said instant n\ 

- the calculation methods before said instant n are restored when time n+AT, 
elapses. 

For what said above, the Cholesky factor can be submitted to triangular 
factorization in this time interval, as follows: 

'L(« + Ar) 0 
U(n) L»_ 

With all the above considerations relevant to the possibility to use in the 
calculations either directly the matrix R(w), or preferably the Cholesky factor L(n) of 
the same, the different steps of the adaptive estimate of the sole vector w^(n) of 
weights characterizing the method embodiment, are connected to the corresponding 
20 ones of the adaptive algorithm as for the utilization method in the computation of low 
delay decided data {*'„}. Said adaptive algorithm in its turn, referred to the detailed 
description of the algorithm developed by the FUNZ block included in the receiver 
sketched in fig. 6, this description continuing to be valid making the modifications 
relevant to the time progress AT of matrix R w (n + AT) or of the relevant Cholesky 

25 factor L(/f + AT) . 

This assumed, the weight vector w opf (w) can be calculated either solving the 
(78) or, preferably, through Cholesky factorization and subsequent solution of the (56) 
which in this new context becomes: 

L*0f+jr)w e » = V'{nfS m (56') 



15 L(«) = 
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Concerning this second opportunity however, it is necessary to analyze how to 
update the Cholesky factor, when it assumes the expression of the (79), applying the 
unitary transformation Q(n) defined in the (76). In the case of use of the embodiment 
contrarily to the conditions that applied to the (61) and (76). it is necessary to update, 
5 at instant n , only the L factor of the spatial covariance matrix R^. Analyzing the 

mechanism according to which the unitary transformations are applied, we can realize 
how it is possible to update matrix L only "delaying" the updating of the other factors. 
It was previously said, making reference to (61) that a unitary rotation is generally 
implemented as a sequence of elementary transformations, Q 0 ,Q l9 - usually Givens 

10 rotations. The updating delay can be introduced remembering what expressed in the 
(62) and (63), where the single elementary transformation Q, acts on the single 
column of factor L and on the column of generator g„. If we apply the 

transformation Q. to the first M elements of the Mh column, where M is the order of 
matrix L, it is possible to delay the updating of factors U and L x , storing the 

15 transformations occurred in L and restarting the updating of U and L JT once low 
delay decisions are newly reliable, that is at the end of the setting period AT of these 
matrices. 

Fig. 19 shows two trends of BER along the burst of a same number of 
demodulated signals, for a strong interferer disturbing up to the 50th symbol. A first 
20 trend (dashed line) is relevant to the out-of-midamble adaptive method of the subject 
equalization method without the characteristic embodiment for rejecting strong 
interferes, while the second trend (solid line) relates to presence of this embodiment. 

E C 

The input signal shows — = IQdB and — = -20 dB ; the trends are averaged over 

N 0 I 

500 bursts. It can be clearly noticed that the new method starts the adaptation first, 
25 thus succeeding in containing the peak in the trend of the error probability, and 
assuring a significant improvement in the quality of the demodulated signal. 

Further bibliographic references useful to the comprehension of the 
"Displacement Theory" applied to the time-variant channel are reported in the following 
papers: 

30 • "Time-Variant Displacement Structure and Triangular Arrays", by: A. H. Sayed, H. 
Lev-Ari and T. Kailath, published on IEEE Trans, on Signal Processing, vol. 42, No. 
5, May 1994. 
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• "Fast Computation of Channel-Estimate Based Equalizers in Packet Data 
Transmission", authors: N. Al-Dhahir, J. M. Cioffi, published on IEEE Trans, on 
Signal Processing, vol. 43, No. 11, November 1995. 

• "A Computationally Efficient FIR MMSE-DFE for CCI-lmpaired Dispersive 
5 Channels", by N. Al-Dhahir t published in IEEE Trans, on Signal Processing, vol. 45, 

No. 1 f January 1997. 

The teachings of the equalization method scope of the present invention, find the 
relevant application in the equalizers included in the functional blocks of figures 6 and 
12. The implementation of the above mentioned blocks is known to those skilled in the 
io art, as well as the implementation of the functions indicated in the text which do not find 
corresponding blocks in the figures, except for the FUNZ block and the relevant internal 
blocks PESI and CAMP. Moreover, as a consequence of the different interrelation 
among the blocks that a joint estimate of parameters w op ,,h opf imposes versus a 

separate estimate, also the equalizer architecture does not result known. 

15 in the detailed description of the FUNZ block below, reference shall be made to 

the sole representation of fig. 12, which enables to implement a sole equalizer capable 
to operate both during the midamble, and out of it, in this last case according to the 
adaptive method, or according to that of the alternative embodiment, after appropriate 
hardware/software selections controlled by events such as the end of the midamble, or 

20 the sudden arrival of a strong interferer. 

As far as the circuit realization of the FUNZ block is concerned, it is preferred to 
have recourse to circuit structures known as "systolic arrays", represented in figures 
21, and 22. These structures have a high parallelism degree, that makes them suitable 
in themselves to operate on the matrices, but mainly to make almost simultaneously a 

25 number of calculations which should otherwise be made in sequence. A 
microprocessor based implementation would be well less effective, just because a 
common DSP is governed by a time sequence of elementary instructions making it 
more suitable to solve serial problems, contrarily to the systolic array that performs all 
the calculations in parallel in one symbol time of 3,69 ^is (or in a double time in case of 

30 joint operation on the semibursts). Even if one wants to have recourse to a DSP of the 
RISC type, the latter would not have a parallelism comparable with that of the systolic 
array, which seems to be therefore "problem oriented", in the practice, the FUNZ block, 
as well as the remaining blocks of the equalizer of fig. 12, are profitably implemented 
through digital integrated circuits of the ASIC type (Application Specific Integrated 
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Circuit). 

Making reference to fig. 20, it can be noticed that the FUNZ estimate block 
includes five differently configured systolic arrays, indicated with SYSTOL1. SYSTOL2, 
SYSTOL3, SYSTOL4, and SYSTOL5, respectively. The vector x{n) referred to the 

5 samples of the training sequence {x n } stored in the receiver reaches a first input of the 
SYSTOL1 block, a second input is reached by vector y(n) , and the third input by a 
MID control signal resulting true only coinciding with the midamble. On a first output 
bus the Cholesky factor T{n) indicated in (61) is present, and reaches a first input of a 
two-input selector SEL, while a second output bus CONT connects SYSTOL1 to a first 

10 input of SYSTOL2. 

The second input of the SYSTOL2 block is reached by vector £'(n) referred to 

the data of the low delay decided sequence {x\} ; a third input is reached by vector 

y(n) coming from the RIT delay block; a fourth input BURST is a bus having reading 

access to a memory in which the receives writes the samples y(n) of the reception 

15 burst; a fifth input is reached by the forgetting factor A, a sixth input is reached by a 
PIC signal, we shall describe below, and finally a seventh input is reached by the MID 
control signal. The Cholesky factor L'(n) is present on the output bus as indicated in 
the (76), that reaches the second input of selector SEL, whose selection input is 
reached by the MID control signal. The bus coming out from selector SEL is connected 

20 to one end of the switch block SWT controlled by a PIC signal having duration AT 
coming from a PICSEL block. This last block receives the MID signal and two other 
signals indicated with PIC and TRANS on three relevant inputs, and on this basis 
generates the PIC signal. The TRANS signal is active for a time AT equal to the 
duration of the transient operation foreseen for the method embodiment; it is a 

25 configurable parameter of the equalizer, like the forgetting factor X . The presence of a 
logic, taking care of the generation and control of all the indicated signals, is obvious. 

The PIC signal comes from a PICRIV block receiving the high frequency receipt 
signal directly from the antenna array ARY. The bus coming out from the other end of 
the SWT block is connected to an input of SYSTOL3, from which the intermediate 

30 vector u„ comes out directed to a first input of SYSTOL4, whose other input is 
connected to the input of SYSTOL3. A vector 1) is present at a first output of 
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SYSTOL4, directed towards an additional input of SYSTOL3. The vector h^in) is 

present at a second output of SYSTOL4, directed towards a first input of SYSTOL5, 
whose second input is connected to the output of the SEL block. At the output of 
SYSTOL5 the vector w opt (n) is present. The two blocks SYSTOL3 and SYSTOL4 are 
5 included in a functional block AUTOVET, whose input coincides with that of SYSTOL3 
and whose output coincides with that of SYSTOL4. 

The architecture of the FUNZ block of fig. 20 enables a full operation of the 
equalizer of fig. 12 both during the midamble and out of the same, according to the 
v method scope of the invention and its alternative embodiment. 

10 During the midamble, the MID signal results true, this condition is used to enable 

Q the array SYSTOL1 and disable SYSTOL2, to control the selector SEL in the selection 

%3 of L(«), and determine the closing of the switch SWT. Consequently, the bus of 

r: SYSTOL1 is extended towards the AUTOVET block extracting the portion relevant to 

ffi matrix L, and uses it to solve the (68) and (67), determining the vector h op ,(n) . More in 

nU 15 particular, the array SYSTOL3 calculates the intermediate vector u n according to the 

q (68) applying the "back-substitution" method, while the SYSTOL4 array calculates 

5 h opX n ) according to the (67), again through "back-substitution". The array SYSTOL5 

T: extracts from the bus coming from SEL the portion relevant to matrices L,andU and 

H uses it, together with vector h opt (n) t to solve the (56) and determine the 

( 

20 vector w„„, {n) , always via "back-substitution". 

op* 

During the adaptive operation, the negate MID signal disables the array 
SYSTOL1 and enables SYSTOL2, and controls the selector SEL to the choice of 
L'(rc) . In lack of a change in the logic condition in the PIC'signal, the switch SWT 
remains closed and the bus of SYSTOL2 is extended towards the AUTOVET and 

25 SYSTOL5 blocks; these last treat the relevant matrix components of the Cholesky 
factor L'(n) like what said for L(n). 

Before the SYSTOL2 array starts operating, the transition from true to negate of 
the MID signal controls the transfer of the content of array SYSTOL1 on the CONT bus 
to be copied in the cells of SYSTOL2. The above mentioned operation enables the 

30 updating of the already existing Cholesky factor L(n) according to the new calculation 
method foreseeing the use of the forgetting factor X. If the above transfer were not 
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foreseen, the array SYSTOL2 would newly estimate the factor L'(«) but with lower 
possibilities versus the previous situation. 

When the PICRIV block detects, at an instant n , the presence of a power peak 
on the antennas, higher than the set threshold, the PICSEL block discriminates 
5 between two possible operation conditions. A first condition is that of strong interferer 
in presence of midamble, in this case the PIC signal results false and the estimate 
during the midamble is not disturbed. The second condition is the strong interferer one 
out of the midamble, in this case the true signal PIC controls the opening of the switch 
SWT, and controls appropriate switches included in the array SYSTOL2 that make a 
10 temporary modification in the structure of the array such to render it suitable to operate 
as foreseen in the method alternative. With true PIC the AUTOVET block does not 
receive anything at input and freezes the value of vector h opt at instant n up to 
elapsing of A7\ The SYSTOL5 block solves the (56), always through 
"back-substitution", using the vector h opt and the matrices LandU of the Cholesky 

15 factor L'(n) derived from the above mentioned structural modifications made to the 

array SYSTOL2, determining in this way the vector w opf (n) . 

When time n + AT elapses, the new condition of the PIC signal controls the 
closing of the switch SWT, but the configuration of the array SYSTOL2 remains as 
modified up to the end of the burst, and the equalizer operation returns to the adaptive 
20 method. 

Concerning the circuitry realization of the different blocks of fig. 20, the following 
must be pointed out: 

- blocks PICRIV, PICSEL, RIT, SEL, and SWT, on the basis of the functional 
description supplied, can be implemented applying the common technical 

25 knowledge; 

- the systolic arrays SYSTOL1, SYSTOL3, SYSTOL4, and SYSTOL5 can be 
represented by an architecture of cells common to that of all known systolic arrays 
of triangular form. Their description therefore requires to specify the nature of the 
cells, to highlight where the array diversify among them and from those of the prior 

30 art. A functional description of the cells indicating the type of arithmetical operations 
made on the entering signal and on the pre-existing content to obtain the new 
content and the output signal, must be considered sufficient to their realization 
availing of the common technical knowledge in the field of systolic arrays and gate 
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arrays. 

- the systolic array SYSTOL2 avail of the same structure of SYST0L1, but includes 
FIFO memories and electronic switches in the descent chain, as we shall see in fig. 
22. The ways to obtain the above mentioned FIFOs and relevant switches in chips, 
5 are still those that enable to obtain the different arrays, that is the gate array 
technology of the ASIC type. 

Of course, the FUNZ block can be implemented through architectures differing 
from the indicated one, it is possible for instance to use a unique systolic array 
grouping the SYSTOLI and SYSTOL2 functions, in this case it should be useless the 
10 transfer of the content to the end of the midamble. This new array must however be 
able to modify at the end of the midamble the arithmetic operations defining the content 
of the cells, since now it must be considered the forgetting factor A . 

With reference to fig. 21, in which the systolic array SYSTOL1 is represented in 
functional blocks, we can notice that there are single blocks having square shape, 
15 interconnected among them, denoted with u i e /,., and peripheral blocks having 
circular shape c t connected to the first ones. The structure of the array has triangular 
shape, as that of the matrix L(n) it reproduces. Actually, each block l n u i , . and 

represents an . element of this matrix. As it can be noticed, blocks l n c i are grouped in 
the area of the upper sharp vertex corresponding to the position of the triangular matrix 
20 L , blocks I s , , c, are grouped in the area of the lower sharp vertex corresponding to 

the position of the triangular matrix L Jf and finally blocks u t are grouped in the 
remaining square area (or, more generally, rectangular) corresponding to the position 
of matrix U . Circular blocks c t form the closing diagonal of the array, and the two 
blocks at the ends indicate the starting and the end of the updating. For ease of 

25 representation 3x3 matrices are shown. 

Fig. 22 indicates the architecture of array SYSTOL2, obtained from that of fig. 21 
interrupting the output line deriving from the blocks corresponding to the elements of 
the last line of the triangular matrix L , including identical FIFO memories having length 
such to introduce a delay equal to the transient AT , and connecting the outputs of the 

30 FIFOs to relevant blocks corresponding to the elements of the first line of matrix U . 
The array SYSTOL2 can operate according to the adaptive method and to the method 
embodiment, the means enabling this possibility shall be described. 
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Fig. 23 shows an electronic switch COM, identical to the other ones included in 
the array SYSTOL2, inserted between two points A and B in fig. 22, corresponding to 
the input and output ports of a generic FIFO. The central point of the switch COM is 
connected to point A of origin of the signal to switch, while the FIFO is connected 

5 between the two output terminals of COM, one of which is connected to point B. The 
active PIC signal controls the switching in order to insert the FIFO between A and B as 
shown by the position of fig. 23; the signal propagating on the corresponding column is 
thus delayed by AT . However, at the end of the transient AT , the FIFO remains 
inserted up to the end of the burst, instant at which occurs the opposite switching that 

10 shortcircuits points A and B, and the delay is excluded returning to the adaptive 
operation in presence of a new burst. In the figure it is also highlighted a clock signal 
CKT at symbol frequency that sets the operation time of the FIFO. There will be as 
many COM blocks as are the columns over the FIFOs. 

Referring to fig. 24, we can. notice that circular blocks are affected by an input 

15 signal in and by an output signal p. Also square blocks are affected by the same 
signals, but they have also an additional input signal p and an additional output signal 
out An exception are the square blocks of the last horizontal row and the last circular 
block of the same, where the output signal p is missing. The content of blocks, 
similarly for each type of block, is denoted with a. 

20 SYSTOL1 and SYSTOL2 arrays are arranged in the figures in such a way that 

input signals, corresponding to the updating vector: 

enter from the left side. More in particular, numeric signals y x (n\... 9 y t 

are directed towards the blocks of the first vertical line of the array, corresponding to 

25 the elements of the first column of the triangular matrix L , while numeric signals 
x, (n),...,x, (*),.. are directed towards the blocks corresponding to the 
elements of the first column of matrix (rectangular) U. Concerning the SYSTOL2 
array, the vector x{n) both during the adaptive operation which during the alternative 
coincides with vector x^n); vector y(n) during the adaptive operation shall be that 

30 coming out from the RIT block, that is y(n - D) ; during the variant on the contrary it 
shall be a vector y(n + AT) consisting of samples 
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y { {ri + AT),...,y ( (n + AT),....,y M (n + AT) of the stored burst, transferred in 
SYSTOL2 and there used in the advanced updating of the Choiesky factor 
L(n + AT) of the spatial covariance matrix R^Cn + Ar) . 

Even if not expressly indicated, provisions are made for a general clock signal of 
the array of figures 21 and 22, at a frequency much higher than the symbol one, which 
simultaneously reaches all the blocks l n w. , 7 Jt . , and c. to control the relevant single 

operation, the bench of FIFO memories of fig. 22 is updated at symbol time. 

The geometric arrangement of the blocks shown in figures 21 and 22 is observed 
on the chip area, where the blocks form a same number of cells; this enables the signal 
originating in the upper vertex to propagate along the descent chain. It is obvious that 
the correct content of the cells shall be available only after a time sufficient to the 
complete propagation of the signal towards the opposite end of the diagonal made of 
the circular blocks, following the longer path between the cells. The single signals 
entering the left side of the array are associated to a corresponding number of arrows 
that indicates the delay according to which their effect develops. 

In the operation, systolic arrays SYSTOL1 and SYSTOL2 are respectively 
referred to the updating of the Choiesky factor L(n) of matrix R(n) as indicated in the 
(61 ) and (76), we repropose for ease of reading: 

[L(„) 0] = [L(n-l) g„]Q(n) (61) 

[L(h) 0] = [VIl(h - 1) g„]Q(«) (76) 
During the strong interferers the updating indicated by the (76) is modified by 
the introduction of the FIFO as seen in figures 22 and 23, but the relations we shall 
now indicate will continue to have value on the elements of matrix L(n + AT) . 

The unitary rotation Q(n) is made associating to each cell a relevant operation 
made on input signals and on the already existing content, to obtain the output signals 
and the new content. In particular, for instance, of the (76): 

- in signals propagate in horizontal direction and correspond to the 

j-.W,...,^") e *,(«),->*,(»),--.*» signals at the input of the 

cells of the first column of the array; for the remaining cells indicated by square 
blocks the in signals correspond to the our signals coming out from the previous 
cells in the horizontal propagation direction of the signal; 

- signals p propagate in vertical direction and are generated at the output of each 
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cell, both square and circular, calculating the following expression: 

pm Tlx m 

where X is the forgetting factor, and a is the new content of the cells calculated as 
follows, concerning the sole circular cells corresponding to the main diagonal of 
matrix L: 



a <_ ^A,\ a \ 2 +\in\ 2 (81) 
the new a being pointed by the arrow; 

- for the remaining square cells the new a is determined jointly to the signal oaf that 
propagates in horizontal direction, in fact, the following matrix relation applies: 

10 \a out] 4— [aVX in] , * = . ^ (82) 

1 J 1 't/HJLp "'J 



where the vector pointed by the arrow indicates the new values of a e out . 

Relations similar to the (80), (81), and (82) are valid for the updating according to 
the (61 ), giving X = L 

Concerning fig. 22, the introduction of FIFO memories in the descent chain of the 
15 systolic array is the simple circuit translation of what said describing the method 
embodiment concerning the updating of the sole Cholesky factor L(w +zir). In fact, 
we said that it is possible to delay the updating of U and L, factors storing the results 
of the elementary transformations Q £ . The delay introduced at instant n is then 
recovered at instant n + AT , that is when the estimator STIM is able to supply the 
20 reliable symbols decided and can therefore restart a normal updating of the Cholesky 

factor L(n). 

As for the functional description of the cells belonging to the systolic arrays 
SYSTOL3, SYSTOL4 and SYSTOL5, that is the mathematical relations developed in 
the cells on the basis of the previous content and of new inputs to generate new 
25 outputs, the sector technician is enabled to obtain the above mentioned systolic arrays, 
applying the steps of the back-substitution method described in the text to solve the 
relevant equations indicated, and availing of the common technical knowledge in the 
field of systolic arrays. 

Profitably, the equalizer belonging to the receiver of fig. 12 (in practice the whole 
30 block structure shown in the figure) is duplicated to process in parallel the two 
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semibursts at two sides of the midamble. A burst reconstruction circuit placed 
downstream the two equalizers puts in the right time order the data of the two 
estimated subsequences, obtaining a unique sequence {x„}of data transmitted 
estimated at maximum likelihood. 



